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^ ' Abstract 

^ ' The Penrose inequality gives a lower bound for the total mass of a spacetime in terms of 

the area of suitable surfaces that represent black holes. Its validity is supported by the cosmic 
■ censorship conjecture and therefore its proof (or disproof) is an important problem in relation 

fT^ ' with gravitational collapse. The Penrose inequality is a very challenging problem in mathematical 

relativity and it has received continuous attention since its formulation by Penrose in the early 
seventies. Important breakthroughs have been made in the last decade or so, with the complete 
resolution of the so-called Riemannian Penrose inequality and a very interesting proposal to 
address the general case by Bray and Khuri. In this paper, the most important results on this 
field will be discussed and the main ideas behind their proofs will be summarized, with the aim 
of presenting what is the status of our present knowledge in this topic. 
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In the early seventies, Penrose [1351 1136j showed that by combining several ingredients of the "estab- 



iy»^ , lishment viewpoint" of gravitational collapse, an inequality of the form 



o: ^^^\liT- (1) 

Q>^ [ V Idtt 



follows, where M is the total mass and A the area of a black hole. Cosmic censorship is one of the 
fundamental ingredients of the argument, and by far the weakest one. Thus, finding a counterexample 
k>( I of ll]) would very likely involve a spacetime for which cosmic censorship fails to hold. In fact, this was 

; I ' Penrose's original motivation to study the inequality. On the other hand, a proof of a suitable version 

^ i of ([1]) would give indirect support to cosmic censorship. Inequalities of this type are collectively 

termed "Penrose inequalities" (sometimes also "isoperimetric inequality for black holes" [62] ) and 
finding suitable versions thereof and trying to prove them has become a major task in mathematical 
relativity. After a first period of heuristic proofs and partial results, important breakthroughs have 
been made in the last ten years. The aim of this review is to try to explain the problem and describe 
the main approaches that have been followed. 

The first observation to be made is the necessity of replacing the area of the black hole in ([T]) by 
the area of a suitable alternative surface. This is because in order to determine whether a spacetime 
is a black hole, detailed knowledge of its global future behaviour is required. On the other hand, 
cosmic censorship is precisely a statement on the global future evolution of a spacetime. In order 
to have an inequality logically independent of cosmic censorship (although motivated by it) the area 
on the black hole must be replaced by the area of a surface which can be located independently of 
the global future behaviour of the spacetime (for instance, directly in terms of the initial data) and 
which is guaranteed (or at least expected) to have less or equal area than the event horizon that may 
eventually form during the evolution. 

The global setup which supports the validity of ^ is well-known and, in rough terms, goes 
as follows. Assume a spacetime {M , g) which is asymptotically flat in the sense of being strongly 
asymptotically predictable, admitting a complete future null infinity ^+ and satisfying J~(,y+) ^ A4 
(see |160j for definitions). The event horizon Ti, is the boundary of J^(J^+) and it is, therefore, a 
null hypersurface at least Lipschitz continuous. Assume, moreover, that the spacetime admits an 



asymptotically flat partial Cauchy surface with total ADM energy Eadm and which intersects TL on 
a cut S. If 7i is a smooth hypersurface then this cut is a smooth embedded surface which has a well- 
defined area 151. For general event horizons, the area \S\ still makes sense provided it is interpreted 
as its 2-dimensional Hausdorff measure (the Hausdorff measurability of S is demonstrated in [37]). 
Consider now any cut Si to the causal future of S along the event horizon. The black hole area law 
[751 [7S] states l^il > \S\ provided the null energy condition holds. This area theorem was proven for 
general event horizons in [37j under much milder asymptotic conditions. From physical principles, 
the spacetime is expected to settle down to some equilibrium configuration. Assuming also that all 
the matter fields are swallowed by the black hole in the process (an external electromagnetic field 
would not alter the conclusions, see Sect. [5]), the uniqueness theorems for stationary black holes 
(see e.g. [85] ) imply that the spacetime must approach the Kerr metric (modulo several technical 
conditions that still remain open, see [IJ for a recent account). For the Kerr metric, the area of 
the event horizon ^Kcr, is independent of the cut (as for any Killing horizon) and takes the value (in 



units G = c = 1) ^Kerr = SttM [KI + ^AP - L^/APj < WttAP where M and L are respectively 

the total mass and total angular momentum of the spacetime (we do not use the more common term 
"J" for the angular momentum to avoid confusion with the energy flux used later). In particular, M 
should be the asymptotic value of the Bondi mass along J^^ . Since gravitational waves carry positive 
energy, the Bondi mass cannot increase to the future [ST] 1139] . Provided the Bondi mass approaches 
the ADM mass Madm of the initial slice (which is only known under additional assumptions, see 
[TT] [74] [Tot] [T68] ). the inequality Madai > V|S'|/167r follows. This inequality is still global in the 
sense that the cut 5* of the event horizon cannot be determined directly in terms of the initial data. 
Penrose's idea was to consider situations in which one could estimate the area of the cut from below 
in terms of the area of some surface that could be located without having to solve the whole future 
evolution. 

One such situation occurs when the initial data set is asymptotically euclidean and contains a 
future trapped surface (see below for definitions). Then, the singularity theorems of Penrose [133] , 
Hawking [7B and others (see [148 ] for a review) state that, provided the strong energy condition holds 
(in some cases the null energy condition suffices), the maximal globally hyperbolic development of this 
data must contain a singularity, i.e. an incomplete inextendible causal geodesic. Very little is known 
in general about the nature of the singularity that forms. More specifically it is not known whether 
the future development admits a complete J^^ and therefore defines a black hole spacetime. The weak 
cosmic censorship conjecture, first proposed by Penrose [134] . asserts that all singularities lie behind 
an event horizon and therefore are invisible to an observer at infinity. Not much is known about 
the general validity of this conjecture, which remains a fundamental open problem in gravitational 
collapse physics, see [161] . Rigorous results are available only in the case of spherical symmetry, where 
the conjecture has been proven for several matter models [36] 150] . 

Under cosmic censorship, the initial data containing a future trapped surface S develops a black 
hole spacetime. A general result on black hole spacetimes [501 HZ] is that no future trapped surface 
can enter into the causal past of J^+. Therefore, the intersection of the event horizon and the initial 
data set S defines a spacelike surface Tis that separates S from the asymptotic region. In general, 
the location of Tis cannot be determined directly from the initial data. Moreover, Tix; can have 
smaller area than S, even though it lies in its exterior. Nevertheless, under suitable restrictions on S 
(details will be given below) it makes sense to consider all surfaces in E enclosing S. The infimum of 
the areas of all such surfaces, denoted by A^in{S), has the obvious property that jTisI > Anun{S). 
Consequently, the inequality 



follows from Penrose's heuristic argument. The need of using the minimum area enclosure of 5' in ([2]) 
was first noticed by Jang and Wald !94] and the first example showing that S may have greater area 
than surfaces enclosing it is due to Horowitz 86J. 

Inequality ^ involves objects defined solely in terms of the local geometry of the initial data 
set and its validity can therefore be addressed independently of whether weak cosmic censorship (or 
any other of the ingredients entering into the argument) holds or not. It is clear that, since >lniin(iS') 
depends on the hypersurface S containing S, we can still take the supremum of the right-hand side with 





(2) 



respect to all asymptotically flat hypersurfaces containing S, and the resulting inequality still follows 
from Penrose's heuristic argument. This gives a clearly stronger inequality. However, it depends on 
the piece of spacetime available and therefore looses the desirable property of depending solely on 
objects defined on the initial data set. 

If the initial data set is asymptotically hyperbolic instead of asymptotically flat (i.e. such that it 
intersects future null infinity in its Cauchy development), the same heuristic argument gives ^ with 
the Bondi mass replacing the ADM mass in the left hand side. 

Another setup where a local geometric inequality is implied by the global heuristic argument of 
gravitational collapse appears in the seminal paper by Penrose |136) and consists of a null shell J\f 
(with compact cross sections) of collapsing dust in Minkowski spacetime. As described in more detail 
below, for any given shape of the shell (restricted to be convex at each instant of time in order to 
avoid shell crossings in the past) and any chosen cross section S on the shell, the energy density of the 
collapsing dust can be adjusted so that is a marginally trapped surface with respect to the geometry 
of the spacetime left after the shell has passed. Since a singularity will definitely form in this setup 
(the shell has self-intersections in its future), cosmic censorship predicts the formation of a black hole. 
Similarly as before, the marginally trapped surface S cannot enter the causal past of J^"*" [45]. The 
intersection of the event horizon with the shell Af must therefore lie in the causal past of S. However, 
since the shell is convex and collapsing, the event horizon cut automatically has at least the same area 
as S. Consequently, the heuristic collapse argument implies M > ^/\S\/l&Tr, where M is the mass 
of the shell. By energy conservation, this can be computed directly on 5* in terms of its geometry as 
a surface in Minkowski. Thus, a geometric inequality is obtained for a class of spacelike surfaces in 
Minkowski spacetime. The status of this version of the Penrose inequality will be discussed in detail 
later on. 

All versions of the Penrose inequality have the structure of a lower bound of the total mass of 
the spacetime in terms of the area of suitably chosen surfaces. Therefore, they can be regarded as 
strengthenings of the positive mass theorem, which says that the total mass of an asymptotically flat 
spacetime cannot be negative, under suitable energy and completeness conditions. The positive mass 
theorem also has a rigidity part, namely that the total mass vanishes only for the Minkowski spacetime. 
The Penrose inequality conjecture also has a rigidity part, which in rough terms states that equality 
will only happen for the Schwarzschild spacetime. Heuristically this can be understood because, in 
the case of equality, the final mass of the spacetime must coincide with the starting one. Therefore, 
no gravitational waves are emitted in the process. This suggests that the whole configuration should 
be stationary. However, the only stationary, vacuum black hole is the Kerr spacetime, and equality 
happens for this metric only if the angular momentum is zero, i.e. if the metric is, in fact, the 
Schwarzschild spacetime. 

Since the original proposal by Penrose this topic has become an active area of research. However, 
the problem has proven to be a difficult one and relatively little progress was made during the first 
decades. The most important contribution in this period is due to Geroch '61], who observed that a 
suitable functional defined on surfaces embedded in a spacelike hypersurface (^,7) is monotonically 
increasing if the curvature scalar of 7 is non-negative and the surfaces are moved outwards at a speed 
which is inversely proportional to the mean curvature of the surface at each point. This is the so- 
called inverse mean curvature flow. This functional (now called Geroch mass) has the property of 
approaching the ADM energy of the hypersurface (provided this is asymptotically euclidean) if the 
surfaces become sufficiently spherical at infinity. Geroch's original idea was to prove the positive mass 
theorem by starting the inverse mean curvature fiow at a point (where the Geroch mass vanishes). 
This idea was then adapted by Jang and Wald 94J, who noticed that the Geroch mass coincides 
exactly with the right-hand side of ([1]) if the starting surface is connected, of spherical topology and 
minimal (i.e. with vanishing mean curvature). If the flow existed globally and the leaves could be 
seen to approach round spheres at infinity, then the monotonicity of the Geroch mass would imply 
the Penrose inequality in the particular case of time-symmetric initial data sets. The flow however, 
generically develops singularities and therefore the argument could not be made rigorous at the time. 
The only case where the Penrose inequality could be proven to hold was in spherical symmetry 
(irrespectively of whether the initial data set was time-symmetric or not) using the so-called Hawking 
mass, which is a generalization of the Geroch mass when the second fundamental form is not zero. 

In the late nineties, however, two important breakthroughs were made. First of all, Huisken and 
Ilmanen [89| were able to prove that Geroch's heuristic derivation could be turned into a rigorous 



proof. This required deep results in geometric analysis and geometric measure theory. These authors 
established therefore the validity of the Penrose inequality for asymptotically euclidean Riemannian 
manifolds of non-negative Ricci curvature and having a boundary consisting of an outermost minimal 
surface (this is now called the Riemannian Penrose inequality). Although this boundary was allowed 
to be disconnected, the Penrose inequality could only be established for the area of any of its connected 
components. Shortly afterwards, Bray 22 was able to prove the Riemannian Penrose inequality in full 
generality (i.e. in terms of the total area of the outermost minimal surface, independently of whether 
this is connected or not). Bray's argument is completely different to the previous one and uses a 
deformation of the given metric in such a way that all the non-trivial geometry gets swallowed inside 
the minimal surface while the total mass of the space does not increase and the area of the horizon 
does not decrease. This process settles down into an equilibrium state given by the Schwarzschild 
metric. Since the Penrose inequality is fulfilled in the final state (in fact, with equality), the Penrose 
inequality holds also for the original space. These two fundamental results have boosted tremendously 
the interest in the Penrose inequality, which has become a very important topic in mathematical 
relativity. Although the general case is still open, several ideas have been proposed to approach it. In 
particular, an important step forward has been made very recently by Bray and Khuri [29] . Although 
several issues remain still open regarding this proposal (in particular, the existence of solutions of 
certain PDE problems) the idea is indeed very promising. 

The aim of this review is to present the most important developments in this field. I have chosen 
not to follow a historical order, and rather I have tried to organize the presentation in a way which I 
consider logically convenient. Since several approaches share some of the techniques, I have collected 
many preliminary results in one section (Sect. While this has the potential disadvantage that, in a 
first reading, it is not clear why and where such results are needed, it has the advantage of simplifying 
the location of the required tools. Readers wishing to enter straight into the topic of the Penrose 
inequality may skip this section and refer back to it whenever necessary. 

An important warning is in order. The topic "Penrose inequalities" is vast and has many ram- 
ifications. Although I will try to cover the most important results in the field, I make no claim of 
exhaustivity. There are several other review papers on this topic in the literature. For the early results 
the reader is advised to look at [119] . The breakthroughs of Huisken & Ilmanen and Bray triggered the 
publication of a number of interesting reviews, where the proofs of the Riemannian Penrose inequality 
were treated in detail, see [87], [88], [83], [24], [25], [26], [145], [124] . 

The structure of the review is as follows. In section [2] the basic facts that will be needed later are 
introduced. This section is divided into five subsections. In subsection 12.11 the geometry of (n — 2)- 
dimensional surfaces as submanifolds of an n-dimensional spacetime is described, and several types of 
surfaces are defined. In subsection 12.21 some relevant variational formulas for the null expansions are 
summarized. In subsection 12. 31 the geometry of codimension-two surfaces as submanifolds of spacelike 
hypersurfaces is discussed and two important existence theorems for outermost surfaces with suitable 
properties are reviewed. In subsection 12.41 the Hawking and Gcroch quasi-local masses, for surfaces 
in four dimensional spacetimes and their general variation formulas are discussed. Subsection 12.51 
recalls the notion of asymptotically euclidean. The following Section (Sect. [3|) is devoted to the 
discussion of the various formulations of the Penrose inequality that have been proposed. Section |4| 
deals with the spherically symmetric case. In Section [5] the so-called Riemannian Penrose inequality 
is treated. The main ideas behind the remarkable proofs of Huisken and Ilmanen and of Bray are 
discussed in subsections 15.31 and 15.41 respectively However, previous approaches based on spinors and 
isoperimetric methods which admit interesting generalizations are also covered in subsections 15.11 and 
15.21 Section [6| discusses the Penrose inequality when the initial data set is asymptotically hyperbolic 
instead of asymptotically euclidean. Section [7| is devoted to describing different attempts that have 
been proposed to address the Penrose inequality in the non time-symmetric case. Particular attention 
is paid to the very promising recent ideas put forward by Bray and Khuri [29| . Section [5| treats several 
strengthenings of the Penrose inequality when particular matter fields or symmetries are present. 
Section [9] is devoted to describing some applications of the Penrose inequality. The paper finishes 
with a some concluding remarks in Section 1101 



2 Basic facts and definitions 



A spacetime {Ai ,g) is a connected Hausdorff manifold endowed with a smooth metric of Lorentzian 

signature (with sign convention {— ,+,H }). Smoothness is assumed for simphcity, many of the 

results below hold under weaker differentiability assumptions. Wc further take A4 oricntablc and 
[M , g) time-orientable and that an orientation for both has been chosen. In this review we will 
be mainly concerned with four-dimensional spacetimes, although higher dimensional results will be 
mentioned at some places. A4 will always be four-dimcnsional unless explicitly stated. 

The Penrose inequality involves the area of codimension-two surfaces. Let us therefore start with 
some basic properties concerning their geometry. 

2.1 Geometry of codimension two surfaces 

S will denote a compact, embedded, oriented, codimension-two surface in an n-dimensional spacetime 
{Ai,g) defined via an embedding ^ : S ^ A4 {S will usually be identified with its image). The induced 
first fundamental form on S*, denoted by h, is assumed to be positive definite, i.e. S is spacelike. Such 
an object will be called simply "surface". The area of S is denoted by \S\. 

At p G 5 we denote by TpS and NpS the tangent and normal spaces of S. This implies TpM = 
TpS (B NpS . For any vector V a.t p this decomposition defines a parallel and a normal vector according 
to V = F'I + V-^. The second fundamental form vector of 5 is defined, as usual, as K{X,Y) = 
— (V^j^y)^ where X and Y are tangent to S and is the covariant derivative on (A4,g). This 
tensor is symmetric and its trace H = trhK defines the mean curvature vector. The trace-free part 
of K will be called 11 in the following. 

S being oriented and the ambient spacetime being oriented and time oriented, it follows easily 
that there exist (globally on S) two linearly independent null vector fields 1+ and l~ orthogonal to 
S. These vectors, which span the normal bundle NS = UpNpS, will always be chosen to be future 
directed and satisfying (Z+ ■ l^) = —2 (dot stands for scalar product with the spacetime metric g). 
They are uniquely defined up to rescalings Fl~^, l~ F~^l~ {F > 0), plus interchange /+ ^ /~. 
The null expansions of S are defined as 6± = {H ■ l^) and they contain the same information as H 
since H = {0-1+ + 6*+/"*-). 

A surface S can be classified according to the causal character of H . Unfortunately, there is no 
unique and generally accepted agreement on how to call the various types of surfaces which arise. In 
this paper, the following notation will be used (in our convention, the zero vector is both future null 
and past null): 

Future trapped: H is timelike and future directed everywhere (equivalently (9+ < and 
e_ <0) 

Weakly future trapped: H is causal and future directed everywhere (^+ < 0, ^_ < 0). 

Marginally future trapped: H is proportional to one of the null normals 1+ with a non- 
negative proportionality factor {6+ = and 9- < 0, or viceversa). 

Past trapped, weakly past trapped and marginally past trapped are defined by reversing all 

inequalities. 

All these types of surfaces have a mean curvature vector of definite causal character and time 
orientation. For brevity, when no indication to future or past is given, then future should be understood 
(i.e. a "trapped surface" is meant to be a future trapped surface, and similarly for the other cases). 

It is often useful to consider surfaces for which one of the null normals can be geometrically selected. 
This preferred normal will be called "outer null normal" and will always be denoted as l+. Notice 
that "outer" here does not necessarily refer to any notion of exterior to S, so this definition must 
be used with special care. Nevertheless, when the null normals can be geometrically distinguished, 
the corresponding null expansions are also geometrically distinct and the following definitions (which 
place no restriction on 6-) become of interest. 

Weakly outer trapped: 6+ < 0. 



Marginally outer trapped or MOTS: 9+ = (equivalently, H points along the outer null 
normal) . 

Weakly outer untrapped: 6*+ > 0. 
Outer untrapped: 6+ > 0. 

These definitions depend not only on the choice of outer direction, but also on the time orientation 
of the spacetime. If the time orientation is reversed (without modifying the outer direction) then 
—l~ becomes the future outer null direction. Therefore by adding the word "past" to any of the four 
definitions above, the surface is meant to satisfy the same inequalities with 9+ replaced by —9-. For 
instance, a past weakly outer trapped surface satisfies 6- > 0. 

Very recently Bray and Khuri [29 have proposed a version of the Penrose inequality which involves 
a type of surfaces whose definition is insensitive to reversals of time direction. However, they still 
require a preferred outer direction. 

Generalized trapped surface: At each point, either < or > (or both). 

Generalized apparent horizon: At each point, either — and < 0, or 0_ = and 

0+ > 0. Notice that generalized apparent horizons are in particular generalized trapped surfaces. 

In the detailed classification of surfaces according to its mean curvature vector presented in [147j . 
generalized apparent horizons are termed null untrapped, which, in my opinion, has the advantage of 
being more descriptive. In fact, my preferred term for these surfaces would be null outer untrapped 
because this emphasizes the fact that they have a privileged "outer" direction. However, notation is 
always a very personal matter and, in this review, I will stick to the name generalized apparent horizon 
put forward by Bray and Khuri in their new proposal of the Penrose inequality. 

In terms of the mean curvature vector, a generalized trapped surface is one where H is allowed 
to point everywhere except along a spacelike outer direction (defined to be a spacelike direction with 
positive scalar product with /+). A generalized apparent horizon has a mean curvature vector which 
is null everywhere and moreover, has non-negative scalar product with any outer spacelike direction. 
Notice that the second condition is ensured provided the scalar product with one outer spacelike 
direction is non-negative. 

Generalized trapped surfaces are indeed generalizations of the previous concepts. In particular, 
weakly future trapped and weakly outer trapped surfaces are automatically generalized trapped sur- 
faces, and the same is true for marginally outer trapped surfaces. 



Name of surface 


Null expansions 


Future trapped 


9+ <0 and e*- < 


Weakly future trapped 


9+ <0 and e*- < 


Marginally future trapped 


[9+ = and 9^ < 0) or {9+ < and 9+ = 0) 


Weakly outer trapped 


9+<0 


Marginally outer trapped (MOTS) 


9+=0 


Weakly outer untrapped 


9+>0 


Outer untrapped 


9+>0 


Generalized trapped surface 


At no point, 9+ >0 and 9- <0 


Generalized apparent horizon 


At each point, 9+9^ = and 9+-9^ >0 



Table 1: Types of surfaces according to their null expansion(s). 



2.2 First order variations of the null expansions 

An important technical tool that will be used often below is the first order change of the null expansions 
9± under variations of the surface. In this section we write down the variations of 9± along the null 
normals l^, which we write as £j-±9±. For the sake of completeness, let us first recall briefly how 



A=0 



variations of a geometric object F are defined. Any vector field ^ in a neighbourhood of S defines 
a variation of the surface as follows. For small enough A e M, let (f>A : S" — > be the embedding 
defined by moving p £ S a parametric amount A along the integral curve of ^ starting at p. We 
write S\ = ^\{S). Assume, for definiteness, that F is a covariant tensor geometrically defined on any 
surface and take a variation $a of S. Define F\{p) = $^(F|$^(p)), i.e. the pull-back of the tensor F 
attached to S\ at the point $a(p)- This defines a curve of tensors at each point p G S. The geometric 
variation is simply £^F = d\F\{p)\\=o. This derivative only depends on the values of ^ on 5* for truly 
geometric objects. However, when F requires additional structure for its definition, then derivatives 
of ^ on S may also arise. This behaviour occurs for instance in £^0+ which requires making a specific 

choice of on each surface 5*^. 

Applying this procedure to the area of S gives (see e.g llOlj ) the first variation of area '^^^^ 

Jg{H ■ ^)ris where ris is the metric volume form of S In particular, the change of area of S for 
variations along the null directions are determined by the integral of the null expansion 6±. 

Let us next write down £j-^9±. Since and Z~ are interchangeable, the corresponding expression 

with also hold for variations along l~ after making the substitution + < — > — everywhere. We start 
with the variation of the induced metric, which is well-known to be 

l+ihAs) = 2/+ • Kab, 
where A, B etc. are tensorial indices on S. This implies 

r+ivs) = e+vs (3) 

for the variation of the volume form on S. For the variation of 0^ along , it is necessary to 
relate the null normal Z+ on each one of the varied surfaces Sx to the corresponding null normal on 
the original surface. As already stressed, this introduces a dependence on the first derivative of Z+ 
via = —1/2{1~ ■ p^l~^)\x=o in the resulting expression, which is the well-known Raychaudhuri 
equation 

£r^ 0+ = Q+e+- K+gK+ - Ric (r+ , r+) , (4) 

where — {l^ ■ Kab) and Ric is the Ricci tensor of {M,g). A more involved calculation gives the 
derivative of 0_ along /+ (see e.g. |2]) 

£j.^e^ ^ -Q+0- - K+gR-"^^ ~ Ric(r+, r_) + iRiem(r+, r_, r+, r_) -l- 2 {DaS^ ^ SaS^) , (5) 

where Riem is the Riemann tensor of (Al, .g), D denotes the Levi-Civita covariant derivative of (S*, K) 
and the one-form Sa is defined by 

5(x) = -i(r_.v^^r+), (6) 

where X is any tangent vector to S. The Gauss identity for S" as a submanifold of (A^, implies 

R{g) = R(K) -{H -11)+ Kab ■ K"^^ - 2Ric(r+, L) -I- iRiem(r+, L, r+, L), 

where R{g) is the scalar curvature of {A4,g). Thus, (O can be rewritten as 

£^^0- = -Q+0- + Km{l+,L) - (i?(/i) -{H-H)^+2 {DaS^ + SaS^) , (7) 

where Ein is the Einstein tensor of g. These expressions are valid in any dimension. The expression 
for £j- follows from ([7|) by interchanging + ^ — and substituting 5"^ — > —Sa (see formula © 
above). 



2.3 Surfaces embedded in a spacelike hypersurface 



Codimension-two surfaces usually arise as codimension-one surfaces embedded in a spacelike hyper- 
surface S of the spacetime M. The induced metric on S will be denoted by jij (Latin, lower case 
indices run from 1 to 3) and the second fundamental form with respect to the unit future normal n 
will be denoted by Aij . The constraint equations relate the geometry of S with some components of 
the Einstein tensor 

i?(7) - A,jA'^ + {ti-^Af = 167rp, (8) 
VjAf - V;tr^A = -SttJi, (9) 

where R{"/) is the curvature scalar of 7, Vi is the covariant derivative in (E, 7), the total energy density 
is defined as Svrp = Ein(ri, n) and the energy flux one-form Ji is defined as SirJiX^ = — Ein(X, n) on 
any vector X tangent to S. The initial data set satisfies the dominant energy condition provided 
P > where \ J\ is the norm of J' with respect to the metric Notice that, despite their names, p 
and J are defined directly in terms of the Einstein tensor. No field equations for are therefore assumed 
either here or elsewhere in this paper. The initial data set is said to be time-symmetrie whenever 
A^j = 0. 

Assume, as before, that S is orientable and that a preferred unit normal rfi tangent to E can be 
selected. In this situation, the null normals will always be uniquely chosen as — n ± rfi. As a 
submanifold of E, S has second fundamental form kab with respect to rfi. Its trace p is the mean 
curvature of 5 in E. The trace- free part of hab will be written as 11™^. 

The decomposition of the second fundamental form Aij into tangential and normal components to 
S will also play a role. We will denote by A^^ the projection of Aij on S, the trace-free part of 
this tensor and q its trace. Then, (see e.g. |101| ) Kab — "^ab" + '^AB'rn, Has = ~n^^n + n™^m 
and H — —qn + prfi, which implies 9± = ±p + q. 

The three remaining independent components of Aij can be encoded in the trace tr^(A) and in 
the normal-tangential components 

Sa = A,jm'e\. (10) 

This definition of 5*^ is consistent with ^ once the choice of described above is made. The general 
variation formulas ^ and ([7]) can be rewritten in this context, after a straightforward calculation 
which uses the constraint equation ([8]), as 

£^p = -lu%U"^^^ - ^p^ - ^Ri^) + lRih), (11) 
£rnq = pA.jm'm^ - A^bk"^^ + DaS'^ + 8Tr{J-m). (12) 

In this 2-1-1 + 1 context, MOTS satisfy p + q = 0, generahzed trapped surfaces satisfy p < \q\ and 
generalized apparent horizons p = \q\. Notice that the latter always have non-negative mean curvature, 
which means that outward variations do not increase the area. This property in fact generalizes to a 
global statement for so-called outermost generalized horizons, as follows. 

Consider an initial data set (E,7y,Ay) with E a compact manifold with boundary and assume 
that the boundary can be split into two disjoint components 9~E and 9"'"E (neither of which is 
necessarily connected). We denote by E° the interior of E, so that E = E° U 9"E U 9+E. We want 
to think of c?~E as the "inner" boundary, which means that we will endow it with a unit normal rn^ 
pointing towards E. Similarly 9"'"E is the "outer" boundary, which we endow with the unit normal 
TO+ pointing outside E (see Figure [1]). 

In many cases, E will contain plenty of weakly outer trapped surfaces. An important question 
arises then: is there an outermost weakly outer trapped surface S7 Intuitively, this means a surface 
which encloses all others, or equivalently, such that no weakly outer trapped surface can penetrate 
outside S. In order to make this precise, it is necessary to have a well-defined notion of "outside" of 
S. This can be achieved by restricting the class of surfaces to those which are homologous to the outer 
boundary. More explicitly, a surface S' is called bounding if it is contained in E° U 9~E (in particular, 
it is disjoint with the outer boundary) and, together with (9"'"E bounds an open domain fl' . For such 



Figure 1: Initial data set with an inner {d S) an outer boundary {d~^T,) and the corresponding choice 
of unit normals. A bounding surface S, and its exterior domain fl, are shown. 



surfaces, the mean curvature p and the outer expansion 9^ is calculated with respect to the normal 
pointing towards H.'. A surface S is outermost in some class if no other surface in the class enters the 
domain 17 bounded by S and The existence of outermost MOTS in this context has been proven 

by Andersson and Metzger assuming the inner boundary to be weakly outer trapped (0+ < 0) and 
the outer boundary to be outer untrapped (0+ > 0), so that they play the role of barriers. The proof 
uses the Gauss-Bonnet theorem so the result is valid in the 3+1 dimensional setting. The precise 
statement is as follows 

Theorem 1 (Andersson &: Metzger [5]) Let (E,7y,yly) by three-dimensional and have an inner 
and outer boundary as described in the previous paragraph. Assume that the inner boundary has 
0-l-(9^S]) < and the outer boundary 9-f-{d^T,) > 0. Then there exists a unique smooth embedded 
bounding MOTS S, i.e. SU — dfl and 0^{S) = 0, which is outermost: any other weakly outer 
trapped surface S' which together with d^T. bounds a domain Q' satisfies Q d fl' . 

An important issue refers to the topology of the outermost MOTS S. A classic result by Hawking 
[80] states that the topology of each connected component of the outermost MOTS has to be either 
toroidal or spherical, provided the initial data set satisfies the dominant energy condition. Recently, 
this result has been extended to higher dimensions by Galloway and Schoen |59j . where it is proven 
that, except for exceptional cases, the outermost MOTS in any spacelike hypersurface of a spacetime 
satisfying the dominant energy condition must by of positive Yamabe type. The exceptional cases 
have been ruled out by Galloway in [60J thus leaving only the positive Yamabe case. In four spacetime 
dimensions, this implies that each connected component of the outermost MOTS must be a sphere. 

An existence result involving generalized trapped surfaces and generalized apparent horizons in- 
stead of weakly trapped surfaces and MOTS has been proven recently by Eichmair ^56j . This result 
does not use the Gauss-Bonnet theorem, so it it not restricted to 3-1-1 dimensions. However, it relies 
on regularity of minimal surfaces, which restricts the dimension of S to be at most seven. Outermost 
generalized apparent horizons have one fundamental advantage over MOTS; they are area outer mini- 
mizing, i.e. they have less or equal area than any other bounding surface fully contained in the closure 
of the exterior region. This area minimizing property makes these surfaces potentially very interesting 
for the Penrose inequality, as first discussed by Bray and Khuri ^29j . This theorem of Eichmair answers 
in the affirmative a conjecture due to H.L. Bray and T. Ilmanen [29] on the existence of outermost 
generalized apparent horizons and its area outer minimizing property. 

Theorem 2 (Eichmair [56j) Let (£,7^^, Ay) be m-dimensional, with 3 < m < 7 and have an inner 
and outer boundary as described above. Assume that the inner boundary is a generalized trapped 
surface p < \q\ and the outer boundary satisfies p > \q\. Then there exists a unique embedded 
generalized apparent horizon S (p = \q\) which is bounding, S U = dU,, and outermost: any 

other generalized apparent horizon S' which together with 9"*"!] bounds a domain fl' satisfies d Q' . 
Moreover S is area outer minimizing with respect to variations in Q. 

In the purely Riemannian case (Aij — 0), MOTS and generalized apparent horizons are simply 
minimal surfaces (i.e. p = 0) and the two theorems above become statements on the existence of an 
outermost minimal surface on a compact domain with barrier boundaries. This particular case was 
already known before (see Sect. 4 of [89! references therein). In fact, the existence of an outermost 
minimal surface played an important role in the proof of the Riemannian Penrose inequality by G. 
Huisken and T. Ilmanen [55] . 



2.4 Hawking mass, Geroch mass and their variation formulas 



Huisken and Ilmanen's proof of the Riemannian Penrose inequality is based on a very interesting result 
by Geroch [BTj who found a certain functional (now called Geroch mass) which is monotonic under 
so-called inverse mean curvature flows (IMCF). This functional is defined for surfaces embedded in a 
Riemannian manifold (E,7). For codimension-two surfaces embedded in a spacetime, the analogous 
to the Geroch mass is the so-called Hawking mass [77] (both coincide when the surface is embedded 
in a time-symmetric initial data set). Since these two masses have similar properties, it is conceivable 
that the Hawking mass may also be useful to tackle the general Penrose inequality (i.e. for non 
time-symmetric initial data sets). A natural first step is to analyze whether the Hawking mass is also 
monotonic under suitable flows. Monotonicity along null directions was first studied by Hay ward |71] . 
More recently, the rate of change of the Hawking mass along an IMCF in an arbitrary initial data set 
(S, '-fij,Aij) was studied in |122j . The null and initial data set approaches can be unified and extended 
by using a spacetime flow formulation, where the two-dimensional surface is varied in spacetime. This 
is discussed in detail in |28j . Let us briefly describe the main points, and how previous results flt into 
this framework. 

For any closed spacelike surface S embedded in a four dimensional spacetime {A4, g), the Hawking 
mass [77] is defined by 



Vs 



(13) 



where x{S) is the Euler characteristic of S. Hence x{S) = 2 for a connected surface with spherical 
topology. For an arbitrary compact surface x('^)=X]2(l — 5) where the sum is over the connected 
components and g is the genus. At this point, the parameter C in p3p is a completely arbitrary 
constant. Its inclusion is relevant for settings where the spacetime has a cosmological constant, or 
when dealing with hyperbolic initial data sets in an asymptotically fiat spacetime (see Sect. |6] below). 
This constant was first introduced in [33j in the time-symmetric context. 

Let V* denote the Hodge dual operation on the normal space NpS (which, under our assumptions, 
has a Lorentzian induced metric and it is orientable, so it admits a canonical volume element once 
an orientation is chosen) . This operation is idempotent and transforms any vector into an orthogonal 
vector with opposite norm. Whenever S admits a notion of "outer" , the orientation of the normal 
bundle will be chosen so that /+* = holds (this implies l^* = —l^)- 

Let us also denote by be the connection on the normal bundle (i.e. if V is normal to S and 
X is tangent to S, V = iy^ x ^)^)- -'^^'^ ^"^Y orthogonal variation vector ^, the derivative of the 
Hawking mass along ^ reads [28] 
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where a is the constant obtained by averaging (^ • H) on S, i.e. 
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is the "transverse part of the 
In particular is non-negative if 



gravitational energy" . This object has good positivity properties [! 
H and ^ are both achronal and have positive inner product. 

From expression (|14p is it straightforward to show that dMff{S\)/dX > whenever (Ein + Cg) 
satisfies the dominant energy condition (recall that a covariant tensor satisfies this property when it 
is non- negative when acting on any combination of causal future directed vectors), H is spacelike and 
the variation vector takes the form 
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where |c| < a is an arbitrary constant. Following Hayward [71] in the null case, these flows have 
been termed uniformly expanding in 28J since the mean curvature along ^ and its dual ^* are both 
constant. 

The case with c — ±a corresponds to a null variation vector and corresponds exactly to the flow 
studied by Hayward [71] . The case c — can be naturally called inverse mean curvature flow vector 



(because ^ 
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H) and monotonicity of the Hawking mass in this case was first mentioned in 



and studied in detail by Frauendiener [S7]. The case |c| < a corresponds to a spacelike flow vector 
and therefore can be rephrased as variations within an initial data set, as follows. 

For a surface S embedded in the initial data set (S, "fij, Aij), the mean curvature vector of S reads 
H = —qfi + prfi and the Hawking mass takes the form 
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For an outward pointing variation vector ^ tangent to S, i.e. ^ 
(|14p can be rewritten in terms of the initial data geometry as 



^TO, the general variation formula 
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where the • operation acts on two tensors of the same class Xai-^-a^i ^Ai - a^ and gives 

{Xa,...a. . = pXa,...a.X^'-''- - qXA,...AY^'-^^- + pYa^-aY^'-^^- . 

This quadratic expression is non-negative provided \q\ < p, or equivalently when 6*+ > 0, 0_ < at each 
point, which is an "untrappedness" condition (note that this condition is basically the complementary 
at each point of the generalized trapped surface condition) . In the particular case of C = and IMCF 
{pe"^ = a) , the variation formula ITBl) was obtained in [122] . 

The first two terms in p6|l are non-negative provided 87r| J| < Snp— C holds (this is automatically 
true if Ein + Cg satisfies the dominant energy condition) and \q\ < p. The last two terms have no 
sign in general. However, the factor in round brackets on each of them integrates to zero. Thus, 
monotonicity can be ensured [15] provided the following conditions hold simultaneously: (i) qe"^ — c 
or is divergence-free and (ii) AsV' — \R{h) + j (jP — g^) = a or e'^p = a, where c and a are 
constants. This leads to four different alternatives for which the Hawking mass is monotonic. The 
uniformly expanding flow condition (|15p with |c| < a corresponds to the case qe'^ = c, and pe'^ = a. 
This follows easily from the expressions H = —qft + pm and its dual H* = —qrn + pn. How restrictive 
are any of these four alternatives on a given spacetime remains an open problem. 

While the Hawking mass is intrinsically a functional on surfaces in spacetime (although it obviously 
admits a rewriting in 3-1-1 language), the Geroch mass is directly a functional on surfaces in a spacelike 
slice (S,7). Its definition is 
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The Geroch mass is fully insensitive to the second fundamental form on the slice (in this sense, it is 
a purely Riemannian object) and always satisfies AIq{S) < Mjj{S). Its variation can be obtained 
formally from expression (|16p simply by putting all terms depending on the second fundamental form 
equal to zero. In particular, this requires substituting Snp — > 1/2R{-^), which follows from the Gauss 
equation for S when the second fundamental form vanishes. Explicitly 
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We emphasize that this result is valid for any spacelike hypersurface of the spacetime, not only for 
time-symmetric ones. For vanishing C this expression was first obtained by Geroch |61| . In the 
time-symmetric case, for arbitrary C and imposing IMCF, this expression was derived in |33j . 



2.5 Asymptotic flatness 

An initial data set (E, 7^^, Aij) is asymptotically euclidean provided T, is the disjoint union of a compact 
set K and a finite union of asymptotic ends , each of which is diffeomorphic to \ i? where B 
is a closed ball. Moreover, in Cartesian coordinates in induced by this diffeomorphism, the 
metric and second fundamental forms have the following asymptotic behaviour (we restrict to 3-1-1 
dimensions for definiteness) 

A., =0(^1), dM,^o{^^, (19) 

where r = ^ SijX'^x^ . The Penrose inequality involves the total mass of a spacetime. In the asymptot- 
ically euclidean setting, the total ADM energy-momentum vector [gj is defined through the coordinate 
expressions 

Eadm ^ lim / {djjij - di^jj)dS\ 
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where Sr is the surface at constant r and c?5* = n^dS with n being the outward unit normal to Sr and 
dS the surface element. This definition depends a priori on the choice of coordinates x^. However, 
they can be shown to define geometric quantities on (Y.,"fij , Aij) provided the constraints satisfy the 
additional decay properties 

i?(7) = o(^), VM; - V,^ = O (1) . (20) 

In fact, the ADM energy- momentum is shown in [T^], generalizing results in [T3], to be well-defined 
under much weaker conditions, where the metric and second fundamental forms belong to appropriate 
weighted Sobolev spaces involving just two derivatives for 7 and one for Aij , and that the left hand 
sides of (l^ni) are integrable on E. The total ADM mass of an asymptotically euclidean initial data set 
is defined as Madm = V^Ida/ - Padm iPadm j- 

The positive mass theorem of Schoen and Yau |1411 1143j states that, provided the initial data 
satisfies the dominant energy condition, p > |J|, then Madm is real and in fact strictly positive 
except whenever (S, jij, Aij) corresponds to a slice of Minkowski spacetime. So far, this theorem has 
been proven in any space-time dimension up to n = 8. Furthermore, it also holds for spin manifolds 
of any dimension, for which Wittcn's spinorial proof [167^ I132j of the positive mass theorem can be 
applied. 



3 Formulations of the Penrose inequality 

The heuristic derivation of the Penrose inequality for a surface S relies on two fundamental facts. 
First and foremost, that the "establishment viewpoint" of gravitational collapse holds and second, 
that S is known to lie behind the event horizon in any situation where a black hole does indeed 
form. This second aspect makes it useful to consider the class of weakly outer trapped surfaces S in 
a given asymptotically euclidean slice. In order to have a useful notion of outer direction in this case, 
one restricts the surfaces to be boundaries of domains I?+ containing the asymptotically euclidean 
end (if there is more than one end, one should be selected from the outset). The outer normal is 
chosen to point inside this domain. Such surfaces will be called weakly outer trapped boundaries 



in the following. Given two surfaces 5*1 and S2 (not necessarily weakly outer trapped) that bound 
respectively exterior domains T)^ and 2?^, we shall say that 5*2 encloses Si provided T)^ contains 
Pj. The complementary of the outer domain T)'^ in E is called "interior region of S"' and denoted 
by V" . It is also useful to define the outer trapped region T-^ in E to be the union of interior 
regions of all weakly outer trapped boundaries. An important fact is that in strongly asymptotically 
predictable spacetimes satisfying the null convergence condition (i.e. Ric(/, /) > for each null vector 
I) the interior region of any weakly outer trapped boundary is contained in the black hole region 
M\ J^[.J^^) (see Proposition 12.2.4 in [160) . and Theorem 6.1 in [45]). Consequently, the same holds 
for 7^. Thus, this type of surfaces are good candidates for a Penrose inequality. As mentioned in 
the introduction, the area of the event horizon cut Tis = n E may be smaller than the area of 
any such surface S. However, for spacetime dimensions n < 8 any surface S that bounds an exterior 
domain, always has a minimal area enclosure, i.e. the outermost of all surfaces which enclose S 
and have less or equal area than any other surface enclosing S. We will denote by ^min(<5') the area of 
the minimal surface enclosure of S. For a weakly outer trapped boundary in a black hole spacetime 
it follows that ylmin(5') < I'^sl, because He encloses S. Hence the heuristic argument by Penrose 
implies Madm > A-cain{S) / IQt^ where S is any weakly outer trapped boundary. Therefore it also 
follows 



AW>supJ%^, (21) 
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where the supremum is taken with respect to all weakly outer trapped boundaries S. This inequality 
can be rewritten in a much simpler way as follows. 

Since a sufhciently large coordinate sphere Sr in the asymptotically euclidean region has positive 
outer expansion, Theorem [T] can be applied |5] to conclude that dT^ is a smooth marginally outer 
trapped boundary, which by construction encloses all other weakly outer trapped boundaries. It 
follows immediately that A-a^i-^{dT^) > Anii„(S') for any weakly outer trapped boundary S, and hence 
also Amin(97^) > sup5 A^i-^[S) . On the other hand dT-^ is itself a weakly outer trapped boundary, 
and hence A,„iii(97^) < sup^ Ai„in(S'). Thus, they necessarily coincide and the Penrose inequality 
(|2ip can be written in the simpler form 



MADM>\h^^i^. (22) 

Although the Penrose inequality, seen as a consequence of cosmic censorship, should be expected to 
hold only for the minimum area enclosure of dT-^ , one can often find in the literature a version of the 
Penrose inequality involving the area of the outermost MOTS, 



This is of course a stronger and simpler looking inequality. However, it is presently known not to 
be true. A counterexample has been found by I. Ben-Dov [2D] by considering a spherically symmetric 
spacetime composed of four regions. The innermost region is a portion of a dust-filled closed FLRW 
spacetime. It is followed by a portion of the Kruskal spacetime of mass M which is then joined 
to another region made of dust, closed FLRW in such a way that the outermost region, which is 
again vacuum, has less mass than the intermediate Kruskal region. In this context, there exists a 
slice for which the outermost MOTS lies in the intermediate Kruskal portion, which leads directly 
to a violation of (|23|) . This example is also a counterexample to a version of the Penrose inequality 
suggested by Penrose himself in ■137], which goes a follows. Suppose that an asymptotically euclidean 
initial data set contains a weakly future trapped boundary S and the spacetime evolves according to 
weak cosmic censorship. Then S is contained in the black hole region that forms. Consider the null 
hypersurface M consisting of past directed null geodesies emanating orthogonally from S towards the 
outer direction, and let k be the corresponding null tangent vector. Since S is weakly future trapped, 
we have 9k > and the area of S does not decrease initially along A/", and this remains to be true 
as long as 9k stays non-negative along the null hypersurface. Since the collapse is taking place in the 



future and M extends towards the past and exterior region, one expects from physical grounds that 
9k does not change sign (this is because M approaches weaker gravitational fields and hence outer 
past directed expansions should be positive). Under these circumstances, the intersection of M with 
the event horizon has at least the same area than 5, and hence the usual heuristic argument implies 
Madai > \/\S\J16tt. The counterexample by Ben-Dov shows that this inequality is not generally true 
either. The difficulty with the argument is double. First of all, the assumption that 6k never becomes 
negative along M need not be true (notice that such a property is not implied by the Raychaudhuri 
equation ([4])). In the example by Ben-Dov, this property is not fulfilled if the starting surface S is any 
of the spherically symmetric future trapped surfaces lying either in the intermediate Kruskal region 
or in the innermost closed dust FLRW region (see Figure 6 in [2D]). Secondly, the argument above 
also relays on the implicit assumption that M intersects the black hole event horizon. In the example 
in |20| these hypersurfaces end up in the white hole singularity and do not intersect the event horizon 
anywhere. 

Weakly outer trapped surfaces on a spacelike hypersurface satisfy 0+ = p+q < and the outermost 
such surface dTj^ is a MOTS (i.e. 6+ =p + q = 0). Given an initial data set (Ej^ij , Aij), a new one 
{Ti,"flj, A[j) can be obtained simply by changing the time orientation, i.e. by setting 7^'^- = jij and 
A[j = —Aij. The heuristic argument to prove the Penrose inequality obviously applies equally well to 
this new initial data set. For surfaces which are boundaries of exterior domains, the future outer null 
direction is now = — L, and therefore the weakly outer trapped surfaces satisfy 6'^ = —9^ < 0, or 
equivalently p — q < (this is, in fact, obvious since the change in time orientation changes q —q)- 
Similarly as before, one can construct the past trapped region as the union of interiors of all such 
boundaries. Andersson and Metzger's result (Theorem [T|) implies again the existence of a unique 
outermost surface dT^^ which satisfies p — q = (i.e. a past marginally outer trapped surface). Thus, 
Penrose's heuristic argument also supports the inequality 



MADM>)J '^"''lfJ^\ (24) 

In a given initial data set, neither 7^ contains 7^ nor viceversa, in general. Although not supported 
by any heuristic argument, a version of the Penrose inequality along the lines of (j23p which has also 
been proposed is 
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Since U 7^ contains both trapped sets, the area of its boundary is not smaller than the minimal 
area enclosures A^i^idT^) or Aynin{dTj7). Thus, ([25]) is in general a stronger inequality than ([22ll 
or ([2^ . So far no counterexample of (|25p has been found and it is not clear whether the inequality 
should hold or not. In the time-symmetric case this inequality reduces to either of (|22|) or (I24p and 
its validity has been proven in general, as discussed in Sect. [5] below. Inequality ((25| also reduces to 
either ([22]) or (|24p in any situation where q is known to be identically zero on the outermost future and 
past MOTS. An example of this behaviour is given by a class of axially symmetric and conformally 
flat initial data sets presented in [52]. For a subclass thereof, the validity of the inequality ([25]) has 
been verified numerically in [97] . 

In the spherically symmetric case, both 7^ and 7^ are obviously spherically symmetric and one of 
them necessarily contains the other. The Penrose inequality ([25]) is known to be true in that case, as 
we discuss in the following section. Since the example by Ben-Dov deals with a spherically symmetric 
spacetime and a spherically symmetric slice, it cannot provide a counterexample to ([25]). On the 
other hand, very few explicit examples are known where (|25p has been verified in situations where it 
does not reduce to either ([22]) or ([24]) . The only cases outside spherical symmetry that I am aware of 
involve a numerical analysis of three special classes of initial data sets |104] . In all these examples the 
inequality ([25]) was numerically confirmed. 

Still another version of the Penrose inequality has been very recently proposed by H. Bray and 
M. Khuri. The fundamental idea behind this proposal is to use generalized trapped surfaces instead 
of weakly outer trapped surfaces. Recall that, on an initial data set, a surface that bounds an 
exterior domain is a generalized trapped surface provided p < \q\. This class of surfaces has two main 



advantages over weakly outer trapped surfaces: first, its definition is insensitive to time reversals, so 
that one can get rid of the complications of dealing with several sets, like 7^, or their union. 
There is only one generalized trapped set Ty, (in a given initial data set) and this contains both 7^ 
and (and, obviously, its union). Moreover, the boundary OTy. is a smooth generalized apparent 
horizon, due to Eichmair's result (Theorem [2] above). The second advantage is that 97s is area outer 
minimizing, i.e. 197s | = ^min(97s). Consequently, a Penrose inequality involving generalized trapped 
surfaces can be stated in a much simpler form in terms of the area of 97s, without the need of invoking 
minimal area enclosures. The Penrose inequality proposed by Bray and Khuri thus reads 



Madm > 




This version is automatically stronger than ([221) or ([M]) because 7s encloses the other trapped sets. 
Arguments in favor of this version will be given below, where we discuss Bray and Khuri's approach in 
more detail. We should emphasize, however, that the heuristic argument by Penrose does not support 
this version. The reason is that one of the key ingredients for the heuristics to go through is that 
the surfaces under consideration lie behind the event horizon provided a black hole does indeed form. 
However, this property is not generally true for generalized trapped surfaces. An example can be 
constructed as follows. Consider a Cauchy slice E of the Kruskal spacetime. The portion of S lying 
in the domain of outer communications (i.e. in the exterior Schwarzschild region) has a boundary 
which necessarily lies in the union of the black hole event horizon, the bifurcation surface and the 
white hole event horizons. Select E such that this boundary has a non-empty intersection with both 
the black hole and the white hole event horizons. Since we are considering the Kruskal spacetime, the 
intersection of E with the black hole event horizon is a smooth MOTS and the intersection with the 
while hole event horizon is a smooth past MOTS, it follows that the outermost generalized apparent 
horizon 97s must contain both of them. However, E can be easily chosen so that the former intersect 
transversally (see Fig[3|). Consequently, some portion of 97s (which must be C^) must lie strictly 
inside the domain of outer communications, i.e. outside the black hole region. 




Figure 2: Spacelike hypersurface E in the Kruskal spacetime which intersects both the black hole 
and white hole event horizons in such a way that the two surfaces defined by these intersections meet 
transversaly. In the figure, only the portion of E lying outside the black hole and the white hole regions 
is shown. However, E is a Cauchy hypersurface, so the intersection with the black hole event horizon 
is a compact surface without boundary, and the same holds for the intersection with the white hole 
event horizon. Since both surfaces are generalized trapped surfaces, the boundary 97s must enter the 
shaded region somewhere. 



4 The Penrose inequality in spherical symmetry 

The Penrose inequality has been proven to hold under the assumption of spherical symmetry, i.e. when 
CEtjjij, Aij) is invariant under an 50(3) action with (generically) S*^ orbits. This inequality was first 
established by Malec and O Murchadha [120j assuming that the initial data is maximal, i.e. ti-^A — 



and in full generality by Hayward [75] ■ In either case the inequality bounds the ADM energy (not the 
ADM mass) in terms of the area of the outermost surface among 5T+ (the outermost MOTS) and 
dT^ (the outermost past MOTS). 

Both arguments are ultimately based on properties of the Misner-Sharp quasi- local energy [129] . 
although this is made fully explicit only in [72] ■ This quasi-local energy is exactly the Hawking mass 
specialized to spherical symmetry. Thus, the inequality in this case can be easily derived from the 
general expressions in subsection 12.41 



For a spherically symmetric asymptotically euclidean initial data set, the outer trapped region 7^, 
being a geometrically defined set, is necessarily spherically symmetric. Its boundary dT^ is therefore 
a metric sphere (if non-empty, which we assume from now on). Now, two different possibilities arise. 
Either the region outside dT-^ contains minimal surfaces (case (i)) or not (case (ii)). In case (i) 
there is an outermost minimal surface S„i lying outside dT-^ . This surface is necessarily a two- 
sphere and therefore encloses dT^ . Let us consider the region outside Sm in case (i) and the region 
outside dT-^ in case (ii). We denote this region by Tiext and the corresponding first and second 
fundamental forms by ^ext and A^xt respectively. This region is free of minimal surfaces, so that 
the metric can be written as ^ext = ^ ImCr) + ^"^ {dO"^ -I- sin^ , with 2m(r) > r. At infinity 

\mir^co^{f) = Eadm- From spherical symmetry, the second fundamental form can be written as 
Aext = W{r)dr'^ + Z{r) {d0^ + sin 
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A direct calculation shows that the mean curvature p and the trace q of the second fundamental 



forms on the surfaces {r — const} are p — ly 1 ~ ^ > 0, g = 2Z/r^. Consequently, the Hawking 
mass ([13)) (with C = 0) of these spheres reads Mnir) = m{r) -\- Z^/(2r). Its radial derivative can be 
obtained immediately from (|16p . using the fact that dr = 1 — H!IiM jfi^ which implies e^p = r/2. 
This gives 

" = 47rr2 U + _jj , (26) 



dr \ T . 

where Jr = (J-dr). The dominant energy condition p > | J| becomes in this case p > Jr\/l — 2m{r)/r 
Furthermore, asymptotically euclidean demands W{r) = 0(l/r^), Z{r) ~ 0{\) so that limr^oo Mnir) 
Eadm- 

Let us first deal with case (i) : since Sm is a minimal surface and cannot be weakly outer trapped (it 
lies outside 97jt), it must have q > and hence it is past weakly outer trapped {6^ > 0). The spheres 
in the asymptotic region have 6^ ^ q — p < 0, and consequently there must exist an outermost sphere 
dTj^ with vanishing 6- (i.e. an outermost past MOTS). In the exterior of dT-^ we have 9+ = p+q > 
and 9^ = —p+q < 0. This implies \q\ < p and thus the bound \Zr~^{l~2m/r)~-^^'^\ < 1 outside dT^ . 
This, together with the dominant energy condition implies that Mnir) is non-decreasing outside dT^ 
(this is just a particular case of the monotonicity properties of Mh discussed above). Being dT^ a 



past MOTS, we have MnidT^ ) ^ i ^ y |/167r. The monotonicity of Mnir) from this surface 
to infinity and the fact that Sm is area outer minimizing establishes the Penrose inequality 



- V 167r - V 167r 

for case (i). For case (ii) we have that dT^ is automatically area outer minimizing and a similar 
argument applies: if there is an outermost past MOTS outside dT^ , apply the monotonicity of 
Mfi from dT-^ to infinity. If there is none, apply monotonicity of AIu from dT-^ to infinity. In either 
case one concludes (since dT-^ is area outer minimizing) , 



Notice that the Penrose inequality does not state (|27p in case (i). There, the minimum area needed 
to enclose dT-^ must be used. This agrees with the discussion in Section [31 As already mentioned 
there, Ben-Dov |20| has found an explicit example in spherical symmetry where the inequality (|27p 



is violated. On the other hand, the argument above in fact proves the Penrose inequahty for the 
outermost of the two surfaces dT^ and dT^ and thus also for 9(7^ U 7^), due to spherical symmetry. 

Notice also that the spherically symmetric Penrose inequality above involves the total energy of 
the slice. This is weaker than the expected Penrose inequality in terms of the total ADM mass. If 
the slice (E,7ij, ^4^) is such that one can generate a piece of spacetime which admits another slice 
with vanishing total momentum, then the Penrose inequality in terms of the ADM mass also follows. 
However, the existence of this piece of spacetime is not always obvious. In any case, it would be 
of interest to find a proof of the Penrose inequality in terms of the total ADM mass in spherical 
symmetry directly in terms of the given data. This might give some new clues on how the general 
Penrose inequality can be addressed. 

5 Riemannian Penrose inequality 

As already mentioned, the field has experienced a fundamental breakthrough in the last decade or so 
with the complete proof of the Penrose inequality in the time symmetric case, first by Huisken and 
Ilmanen |89j for a connected horizon and then by Bray [22j for an arbitrary horizon. Both papers 
dealt with the four-dimensional case. However, while Huisken and Ilmanen's proof is very specific to 
four dimensions (because of the use of the Geroch mass), Bray's approach can be generalized to any 
spacetime dimension not bigger than eight, as recently shown by Bray and Lee |27|. 

By definition, an initial data set is called time-symmetric whenever Aij = 0. This has two imme- 
diate consequences, namely that the ADM three- momentum vanishes identically, so that Madm ~ 
EadMi smd that only one constraint equation remains. 



which gives i?(7) > provided the dominant energy condition holds. In fact, the weak energy condition 
(defined as Ein(M, u) > for all causal vectors) suffices in this case. 

Another immediate consequence is that 0+ = — 0_ = p and, hence, the trapped region 7^, the 
past trapped region T-^ and the generalized trapped region Tj] all coincide in this case. Its boundary 
Sm is the outermost minimal surface, which is non-empty as soon as there is a bounding surface with 
negative mean curvature (with respect to the normal pointing into the chosen asymptotically euclidean 
end). This is a corollary of Theorems [Hand [2] above. However, in this context this result is known to 
hold even in more generality. As discussed in 89J, following classic results on minimal surfaces jl27j . 
it is sufficient to define the trapped region K as the image of all immersed minimal surfaces in E 
together with the bounded connected components of its complementary. It follows that the boundary 
of this set is a collection of smooth, embedded, minimal 2-spheres and that any connected component 
of E \ X is an "exterior" region, i.e. an asymptotically euclidean manifold free of minimal surfaces 
(even immersed) and with compact and minimal boundary composed of a finite union of 2-spheres. 
The Penrose inequality therefore becomes an inequality relating the total mass and the area of the 
outermost minimal surface Sm in (E, 7) with respect to the chosen asymptotically euclidean end. The 
corresponding inequality 



is usually termed "Riemannian Penrose inequality" since it involves directly Riemannian manifolds, 
with no further structure coming from the ambient Lorentzian manifold. 

As already mentioned, the Penrose inequality has a rigidity part, namely that equality is achieved 
only for slices of the Kruskal extension of the Schwarzschild metric. The time-symmetric slices of this 
spacetime define a manifold (Esd, = \ 0,7sch) with induced metric 



The surface r = m/2 is minimal and separates the manifold into two isometric pieces (corresponding 
to the two asymptotic ends of the Kruskal metric). The rigidity part of the Riemannian Penrose 
inequality states that if equality is achieved in (^5)1 then the region in (E,7) outside its outermost 
minimal surface Sm is isometric to the domain r > m/2 of (Es^h,7sch) with m = Madm- 




(28) 




(29) 



Before discussing the breakthroughs of Huisken & Ilmanen and Bray, let us discuss some previous 
attempts to address the Riemannian Penrose inequahty. 



5.1 Spinor methods 

Shortly after Schoen and Yau proved the positive mass theorem, Witten [167j proposed a completely 
different method using spinors (see [132j for a rigorous version of Witten's ideas). It is natural to 
ask whether spinorial techniques can be also applied to the Penrose inequality. After all, the spinor 
techniques had been successfully extended to prove the positive mass theorem in the presence of black 
holes (more precisely, marginally trapped surfaces) in the initial data [51] (see [HI] for a rigorous 
proof). The main difficulty lies in finding suitable boundary conditions on Witten's equation on 
the boundary of the black holes so that the boundary term arising by integrating the Schrodinger- 
Lichnerowicz [146[|111] identity can be related to the area of the black hole. An interesting attempt to 
achieve this in the Riemannian case is due to M. Herzlich ^81], who obtained a Penrose- like inequality 
involving not just the total mass and area of the minimal surface but also a Sobolev type constant 
of the manifold. Depending on the space under consideration, Herzlich's inequality may turn out to 
be stronger or weaker than the Penrose inequality (see |121| for a limiting case where this inequality 
reduces simply to a positive mass statement). Nevertheless, the inequality is still optimal in the sense 
that equality is achieved only for the Schwarzschild manifold. 

The class of manifolds considered in [ST consists of asymptotically euclidean 3-dimensional, ori- 
entable Riemannian manifolds (S,7) having an inner boundary dT, which is topologically an and 
geometrically a minimal surface. No assumption is made on whether this surface is the outermost 
minimal surface in (S,7) or not. This already indicates that the inequality to be proven cannot be 
the standard Penrose inequality because minimal surfaces with large area can be shielded from an 
asymptotically euclidean region with small mass by an outermost minimal surface with sufficiently 
small area, while maintaining non-negative Ricci curvature everywhere. This is often called "shielding 
effect" and explicit examples are easily constructed by cutting the Schwarzschild manifold (|29)l on a 
sphere at r < m/2 and attaching a piece of of large radius with a cap near the north pole re- 
moved (this manifold belongs to the class of initial data leading to the so-called Oppcnheimcr-Snyder 
spherical dust collapse [131 ) 

The inequality proven by Herzlich reads [81j 



where denotes, as usual, the collection of smooth functions with compact support. Equality in 
(|30)l occurs if and only if (E, 7) is isometric to the exterior of the Schwarzschild manifold (|29p outside 
the minimal surface r = m/2. 

The proof is based on an improvement of the positive mass theorem (also proven in [ST]) valid 
for asymptotically euclidean Riemannian manifolds (5],g) of non- negative scalar curvature having an 
inner boundary dH which is topologically a sphere and which may have positive, but not too large, 
mean curvature (with respect to the normal pointing towards infinity). More precisely, the mean 
curvature p must satisfy the upper bound 



This positive mass theorem also has a rigidity part that states that equality is achieved if and only if 
(S, g) is the exterior of a ball in Euclidean space (this implies, in particular, that the only surface S in 
R'^ which is topologically an S*^ and which satisfies p < 4^^/^^/\S\ is in fact a sphere). The proof of this 
theorem involves finding appropriate boundary conditions for the Witten spinor on dT, so that the 
boundary term in the Schrodinger-Lichnerowicz identity on the inner boundary gives a non-negative 
contribution. 




(30) 



where cr is a geometric, scale invariant, quantity on (^,7) defined as 





(31) 



The Penrose-like inequality (|30p is then proven by finding a conformaUy rescaled metric g = f*g, 
with / ^ 1 at infinity in such a way that g has vanishing scalar curvature and the mean curvature 
of 91] with respect to the metric g saturates the inequality (151]) . Moreover, the conformal rescaling 
is shown to decrease the mass at least by the amount given in the right-hand side of (|30p . Here is 
where the quantity a arises. Since the mass after the rescaling is non-negative due to the positive 
mass theorem above, the Penrose-like inequality follows. The rigidity part holds because, in the case 
of equality in (j30|) . (^,7) is conformal to the flat metric outside a ball and, moreover, its curvature 
scalar vanishes (if this was non-zero, then the decrease in mass due to the conformal rescaling would 
be larger than the right-hand side of (j30p , which cannot occur in the equality case) . Since the only 
conformaUy flat, scalar fiat and asymptotically euclidean Riemannian manifold with a minimal surface 
is the Schwarzschild space (f^ . the rigidity part follows. 

This Penrose-like inequality has been generalized in three different ways. First, Herzlich [M] 
extended the results to spin manifolds of arbitrary dimension n. In this case the boundary dT, is 
assumed to be a compact and connected minimal surface of positive Yamabe type (i.e. such that 

it admits a metric of positive constant scalar curvature Rq)- Denoting by y the Yamabe constant 

2 

(i.e. y — _Ro|9I]| where \dT,\ is the {n— l)-dimensional volume of the boundary), the Penrose-like 
inequality reads 



The idea of the proof is similar to the three-dimensional case. 

The second generalization |118j involves maximalimtial data sets (S, 7^ , Aij). Although this result 
is not a Riemannian Penrose inequality, the methods used are very similar to the previous ones and it 
is therefore natural to include it here. The idea of the proof involves, again, a positive mass theorem 
for manifolds with suitable boundary and finding an appropriate conformal factor which transforms 
the data so that the previous mass theorem can be applied while decreasing the mass by a certain 
amount. 

More precisely, the class of manifolds (S, jij , Aij ) under consideration are asymptotically euclidean 
spin manifolds of arbitrary dimension n satisfying the dominant energy condition p> \ J\. As before, 
the boundary is connected, compact and of positive Yamabe type. The condition of being minimal 
is replaced by three conditions, namely (i) the mean curvature p (with respect to the direction pointing 
towards infinity) is non-positive (ii) q > |p| (or alternatively —q > \p\) everywhere, and 



{Sa is defined in(|10p). Under these circumstances, the total energy satisfies the bound [118] 



The constant Q is always non-negative under assumptions (i)-(ii). Thus, inequality ((33)) is weaker 
than the corresponding one in the time-symmetric case ([5^ . The conditions on the boundary 9S 
are somewhat surprising. Conditions (i) and (ii) state that the mean curvature vector of the surface 
points inward (in the sense that its product with the outer normal m is non-positive) and is causal past 
(future) everywhere. Thus, the boundary is indeed weakly past (future) trapped. However, these type 
of surfaces are necessarily not area outer minimizing (except in the very special case that their mean 
curvature vector vanishes identically). Consequently, the minimal area enclosure of the boundary lies, 
in general, inside E. Thus, this Penrose-like inequality is obtained for a surface for which the original 
Penrose inequality is not expected to hold. It is therefore interesting that such an inequality exists. 

The third generalization is due to M. Khuri |108j and again involves a non-vanishing second 
fundamental form. Since the method uses the Jang equation in a fundamental way, we postpone its 
discussion to the end of Subsection 17.41 




(32) 



where a is the analogous scale invariant quantity in higher dimensions 






(33) 



5.2 Isoperimetric surfaces 



An interesting attempt to proof the Riemannian Penrose inequality is discussed in H.L. Bray's Ph.D 
thesis [Sff • The idea is to consider a special class of surfaces which interpolate between the outermost 
minimal surface S"™, tends to large round spheres at infinity and for which the Geroch mass is non- 
increasing. Although this sounds familiar with the inverse curvature flow argument of Geroch, the idea 
exploited by Bray is in fact very different. Indeed, instead of using flows of surfaces. Bray considers, 
for any given volume, an area minimization problem (i.e. an isoperimetric problem). Assume that the 
outermost minimal surface Sm is connected and consider the class of surfaces in the same homology 
class as Sm- One can associate to each surface S in this class the volume bounded between Sm and S 
(counted negatively in the portion where S lies inside Sm and positively where it lies outside). For a 
given value of ^ > 0, consider the collection of surfaces which bound, together with Sm, precisely a 
volume V and define A(V) as the infimum of the corresponding areas. If the infimum is attained on a 
surface Sy, then this surface is obviously of constant mean curvature p{V) because it is the solution 
of an isoperimetric problem. Bray's idea is to show that the Geroch mass is a non-increasing function 

of y. 

Let us start by assuming that Sy exists and that the function A{V) is twice differentiable. Then, 
the first variation of area together with the fact that the area is the first variation of volume, gives 
A'{V) — p{V) (prime denotes derivative with respect to V) and the Geroch mass (fT7|) becomes simply 



provided the surface Sy is connected and of spherical topology. Consider now a variation of Sy along 
its outer unit normal m with unit speed. If we denote by Sy{t) the corresponding fiow of surfaces 
(with Sy{t = 0) = Sy), the second variation of area gives 



d'\Sy{t)\ 



^ {£,np + p') vs. = 1^ (^-in-^n-^^ + _ l^(^) + 



where we have used (|Tl1) in the second equality. Using now the Gauss-Bonnet theorem and the non- 
negativity of the curvature scalar of (S,7) we conclude ^-^^J^\t=o < 4:Tr + \p{V)^ A{V) . Denoting 
by V{t) the volume bounded by Sy{t), it follows that A{V{t)) < 15^(^)1 because A{V) is the infimum 
of all areas bounding a volume V. Since these two functions touch at t = 0, it follows ^it^^*"*"* \t=o ^ 

d^\S (t)\ 

— It^o- Performing a change of variables t ^ V , one concludes 

A"(V) < =^ Mr(V) > 

Although the argument just described assumes that A(V) is differentiable, the conclusion is still valid 
if a suitable weak (distributional) derivative is taken [30]. The other two assumptions that enter into 
the argument are (i) that Sy is connected and of spherical topology and (ii) that the infimum of 
A{V) is in fact attained i.e. that the surface Sy exists. The first condition turns out to be crucial 
and needs to be imposed as an assumption (in fact, the hypothesis can be relaxed to demand just 
that for each > if one or more minimizers of A{V) exists then at least one of them is connected; 
this is termed "Condition 1" in [30]). Most of the technical work in [30] consists in showing that 
condition (ii) (i.e. the existence of a minimizer for all > 0) imposes no extra restriction. To that 
aim. Bray first argues that the class of metrics (^,7) can be chosen to be exactly Schwarzschild at 
infinity (a similar, but weaker, reduction is also used in Bray's full proof of the Riemannian Penrose 
inequality summarized in subsection 15.41 and will be discussed in more detail there). The heart of 
the proof of existence of the minimizer Sy consists in proving that the isoperimetric surfaces in the 
Schwarzschild spacetime are given by the spherical orbits of the 50(3) isometry group. Contrarily 
to what one may expect, proving this fact is not a trivial matter. The fundamental idea behind the 
construction is to use a comparison metric (in this case the flat metric on R^) for which one knows 
that the spheres are the solutions of the isometric problem. This argument, which was tailored for the 
Schwarzschild metric, has been simplified and substantially extended in |23| . where simple conditions 
are found on a given spherically symmetric metric which ensure that the spheres are the minimizers 



of the isoperimetric problem. Summarizing, Bray proves in [30 that the Penrose inequahty (|28p holds 
for each asymptotically euclidean Riemannian manifold which has a connected outermost minimal 
surface Sm and satisfies "Condition 1" . 

5.3 Huisken and Ilmanen's proof 

The heuristic idea behind the proof of Huisken and Ilmanen was first proposed by Geroch [61j and 
is based on the observation that the Geroch mass (fT7|) (with C = 0) is monotonically increasing if 
the surfaces are moved by inverse mean curvature (i.e. pe^ = 1) provided the scalar curvature of 
(S, 7) is non-negative. This fact is clear from (fT8|) since pe^ = a = 1 and the derivative of Mg{S\) 
is then a sum of non-negative terms. Another immediate property of Mq is that its value on any 
connected, topologically 5*^ minimal surface is Mq = ^/\S\J(T6^^) . For surfaces Sr — r — const in 
the asymptotically euclidean end S]°°, the asymptotic decay of 7 implies limr^oo MciSr) = AIadm- 
Thus, Mg indeed interpolates between the left and right-hand sides of the Penrose inequality. 

Geroch original idea was to prove the positive mass theorem by starting the inverse mean curvature 
flow from a point so that Mq — initially (the point can be approached as the limit of very small 
coordinated spheres). The positivity of mass would follow provided the inverse mean curvature flow 
remained smooth all the way to infinity and the flow approached large coordinate spheres in the 
asymptotically euclidean end. Jang and Wald [94| realized later on that the same argument could 
be used to prove the Penrose inequality by starting from the outermost minimal surface. In fact, 
p vanishes on this minimal surface, and hence the velocity e"^ diverges there. However, by starting 
the flow from surfaces which approximate the outermost minimal surface from outside, the Penrose 
inequality would follow. However, it was immediately realized that the flow will not remain smooth 
in general and that singularities will develop. The flrst variation of area implies that, as long as the 
flow remains smooth, the area of the surfaces increases exponentially. Consider now, as a very simple 
example, two disjoint balls in Euclidean space. By symmetry, each surface will evolve under IMCF 
to a larger sphere with exponentially increasing radius. Thus, the two surfaces will necessarily touch 
in finite time and the flow cannot remain smooth forever. A less trivial example, discussed in [87], 
consists of a thin torus in Euclidean space. The differential equation satisfied by p under IMCF is of 
parabolic type and the maximum principle implies that p stays bounded above in terms of its initial 
value (and the background geometry) as long as the fiow remains smooth. Thus, the torus will flow 
outwards at positive speed bounded away from zero and it will thicken. However, sufficiently thick 
torus have vanishing mean curvature at points on their inner rim. The speed becomes infinity there 
and the flows necessarily stops being smooth. 

The presence of singularities made this idea dormant for decades. The only case were the method 
was made rigorous involved a particular case of metrics called quasi- spherical metrics with divergence- 
free shear [17] , where the flow was seen to remain smooth all the way to infinity. Huisken and Ilmanen's 
fundamental contribution was to define the fiow in a suitably weak sense so that the singularities could 
be treated (and, in fact, basically avoided). 

An important ingredient in Huisken and Ilmanen's approach is the use of a level set formulation 
for the IMCF (which is a geometric parabolic fiow). This means describing the leaves of the fiow as 
the level sets of a real function u on E. The IMCF condition translates directly into the following 
degenerate elliptic equation for u 



where diYjV is the divergence of V in (£,7). In principle it is possible that u remains constant on 
open sets, which has the immediate consequence that the fiow may jump across regions with positive 
measure. The fundamental idea is to use these jumps precisely to avoid the singularities that the 
smooth IMCF would otherwise have. In order to achieve this, Huisken and Ilmanen find a variational 
formulation for ((35|) . This equation is not the Euler-Lagrange equation of any functional. However, 
by freezing the right-hand side to |Vm|, the authors write down the functional (which now depends 
on u) 



The critical points of this functional with respect to compactly supported variations of v gives (1351) 




(35) 




(36) 



with u replaced by v on the left-hand side. One then looks for functions u which minimize their own 
functional, thus giving ([35]) . 

This variational formulation has a geometric counterpart which, in rough terms, implies that each 
of the level sets of u (defined as d{u < t} for positive t) is area outer minimizing. Thus, if we 
start from a smooth surface 5*0 which is area outer minimizing and with positive mean curvature, the 
IMCF (which enjoys short time existence thanks to its geometric parabolic character) will evolve the 
surface smoothly for some "time" A. For small enough values of A, each level set will be outer area 
minimizing. However, there may well exist a value for which S\ ceases to be area outer minimizing. 
Take the smallest of such values, Ai. This means that there exists a surface S'^^ which encloses S\-^ 
and has less or equal area. In fact, it must have equal area because if if had strictly less area, it 
would also have less area than some close enough previous leaf (the flow is smooth up to Ai and the 
area changes smoothly). The surface may have common points with S\j^. However, the surface 
S"^^ \S'ai must be a minimal surface because otherwise it would admit a compactly supported variation 
that would reduce its area while still remaining outside Sx-^ and enclosing it. However, this varied 
surface would (i) enclose all surfaces in the flow from 5*0 up to Sx-^ , (ii) have larger area than 5*0 
(because the latter is area outer minimizing) and (iii) have less area that Sx^ ■ So, for some value of 
the flow parameter smaller than Ai it would have the same area, and it would enclose it. This would 
contradict the definition of Ai as being the smallest value with such property. Notice that, a priori, 
this outer minimizing prescription does not exclude that the flow may become singular already in the 
interval < A < Ai, where all the leaves remain area outer minimizing. That this cannot happen is 
one of the several statements that Huisken and Ilmanen had to show, and which makes their proof 
technically difficult. 

In the torus example before, the smooth flow would thicken the torus (in an axially symmetric 
way) until a surface is obtained which has exactly the same area than the surface obtained by closing 
its hole by two horizontal planes. It is clear in this example that the new pieces have vanishing mean 
curvature. 

Since the level set function u is such that each level set is outer area minimizing, this means that 
the IMCF evolves smoothly for as long as the leaves remains area outer minimizing until a surface 
which is not area outer minimizing is eventually reached. At this point, the surface jumps (meaning 
that u remains constant between the two surfaces). The flow should be continued from the new surface 
outwards. Since the new surface has pieces with p = 0, the IMCF cannot be defined there in a classical 
sense. However, the weak formulation in terms of level sets also takes care of this. 

A crucial condition for the monotonicity argument to go through in this formulation is that the 
Geroch mass does not increase in the jumps. We know that the total area does not change. Moreover, 
the new pieces of the surface have p = 0, while the deleted pieces have p > because the flow was 
smooth up to and including Sx- Thus, across the jump the Geroch mass is non-decreasing provided 
the Euler characteristic of the surface does not decrease. The example with the torus shows that in 
general the topology of the surface may change across the jump, in that case changing from a surface 
with genus one to a topological sphere. In this particular case, the Euler characteristic increases 
and the Geroch mass would be monotonically increasing. The question is therefore whether this is 
a general property or not. The example with two spheres in Euclidean space shows that this cannot 
be true in general. There, the value of the Euler characteristic is initially four and after the jump a 
topological sphere forms, for which x = 2. A general result by Huisken and Ilmanen [89] is that a 
topologically S'^ surface in an asymptotically euclidean three-manifold (S,7) outside the outermost 
minimal surface cannot jump, via the weak formulation of the IMCF, to another surface with smaller 
Euler characteristic, i.e. no holes can appear on the surface after the jump. Nevertheless, the outermost 
minimal surface S in (S,7) needs not be connected. Consequently, for the Geroch mass to remain 
monotonic under the weak IMCF, a connected component Si of the outermost minimal surface must 
be chosen as initial surface. The Penrose inequality proven by Huisken and Ilmanen is, therefore. 



where i runs over the connected components of the outermost minimal surface. 

As mentioned before, the Riemannian Penrose inequality also has a rigidity part. In this setting, 
this states that equality is achieved in (|57|) if and only if the region in (S,7) outside its outermost 
minimal surface Sm is isometric to the domain r > m/2 of (Ssch,7sch) (|29p with m — Madm- 




(37) 



Heuristically, it is clear that some rigidity is to be expected already from the monotonicity property 
of the Geroch mass. If equality is achieved, then the derivative of the Geroch mass is everywhere 
zero. In particular, in the smooth part of the flow, -R(7) = 0, tp = const and each leaf is totally 
umbilical (i.e. 11™^ = 0) at each point. The IMCF condition gives p = const, on each leaf. Since 
d\ = e^rh = ^rh, the variation formula pTjl implies ^ = ~jP+ "^l^- This implies that each leaf is a 
metric sphere (because R{h) — const.). The exponential grow of the area gives \S\\ — |S'o|e'''. Defining 
a new variable f = y/\Sx\/AT:, the metric 7 takes the form 7 = ■^^df'^ + f^dfl. The expression for 

the Geroch mass and the fact that it takes a value m independent of A gives f^p^ — 4(1 ~ and 
hence 7 is the Schwarzschild metric outside the minimal surface f ~ 2m. This is, of course, a heuristic 
derivation. Huisken and Ilmanen are able show the same results using only the weak flow. 

Huisken and Ilmanen's proof has interesting side consequences. First of all, the argument does not 
use the positive mass theorem anywhere and therefore the method gives an independent proof of this 
important result. In the presence of an outermost minimal surface, positivity of Madm is an obvious 
consequence of the theorem. If E is free of minimal surfaces, it suffices to start the weak inverse mean 
curvature flow at one point, making rigorous the original idea of Geroch. 

Another interesting consequence deals with the so-called Bartnik mass (or capacity) [15] , defined as 
follows. A 3-dimensional Riemannian manifold (fi,7) with no boundary and compact metric closure 
is called admissible if it can be isometrically embedded into a complete and connected asymptotically 
euclidean three-dimensional manifold (S, 7) satisfying the following three properties: (a) the curvature 
scalar is non-negative, (b) the boundary of E is either empty or a minimal surface, and (c) S is free 
of any other minimal surfaces (the original definition excludes minimal surfaces even at the boundary 
of E, this modification is due to Huisken and Ilmanen [89]). (^,7) is called an admissible extension. 
The Bartnik capacity is the infimum of the ADM masses of all possible extensions of (^^,7). By the 
positive mass theorem, it is immediately non-negative and it was conjectured in |161 1171 118] to be 
positive if ($^,7) is not a subset of Euclidean space. Huisken and Ilmanen used the weak IMCF to 
reach a slightly weaker conclusion, namely that if the Bartnik capacity is zero, then (^^,7) is locally 
isometric to Euclidean space. The idea of the proof is to take a point p m Q where the metric is 
non-flat and choose any admissible extension (E, 7). Consider the weak IMCF starting at p. For some 
small enough value of the flow parameter the corresponding leaf S\g can be proven to have a positive 
Geroch mass, and to lie within Q independently of the extension (this means in particular that S\g 
is area outer minimizing independently of the extension). Thus, monotonicity of the IMCF gives the 
lower bound Madm > ^-tG(S\„). Since this is independent of the extension, the result is established. 

Another consequence of Huisken and Ilmanen's construction is the so-called "exhaustion" property 
of the Bartnik mass. This can be formulated as follows. Consider an asymptotically euclidean Rie- 
mannian manifold (E,7) and any increasing collection of bounded sets Sli C E satisfying Uifli — E. 
Then the Bartnik mass 1713(^1) tends to the ADM mass when i ^ 00. This is proven in "Mj by taking 
suitable coordinate balls within each fi^ in such a way that the coordinate radius Ri tends to 00. It is 
shown, by using the IMCF, that the Bartnik mass of each fli is bounded below by the Geroch mass 
of the corresponding coordinate balls. The Geroch mass of the coordinate balls tends to the ADM 
mass when i?i — > 00. Since the Bartnik mass of Qi is obviously bounded above by Madm of E, the 
exhaustion property ms(f^i) Madm follows. 

5.4 Bray's proof 

H.L. Bray was able [22 to prove the Riemannian Penrose inequality for arbitrary outermost minimal 
surfaces, with no restriction of connectedness. Bray's proof also uses geometric flows in an essential 
way, but instead of flowing surfaces in a fixed Riemannian manifold, as in Huisken and Ilmanen's 
proof, it uses a fiow of metrics. The idea is to modify the metric with a one-parameter family of 
conformal factors so that the curvature scalar remains non-negative, the new horizon with respect to 
the conformal metric has the same area as the original one and the total ADM mass does not increase 
during the process. If, moreover, the fiow can be defined in such a way that, as i — > 00, the Riemannian 
manifold outside its horizon converges to the r > to/2 portion of the Schwarzschild manifold (j29p in 
a suitable sense, then the Penrose inequality for the original manifold follows because the inequality 
(in fact equality) holds in the limit. This conformal fiow also allows Bray to show the rigidity part of 
the Penrose conjecture, namely that equality holds if and only if the starting Riemannian manifold is 



isometric to the exterior region of Schwarzschild. 

Bray's statement of the Penrose inequahty deals with three-dimensional asymptotically euclidean, 
Riemannian manifolds (£,7) with one or several asymptotic ends (in the latter case, one of them is 
selected and the mass and the concept of "outer" are taken with respect to this end) and which contain 
an area outer minimizing horizon 5*0. By definition, a horizon is a smooth, compact, minimal surface 
which is a boundary of an open set f2~ containing all possible asymptotic ends on (£,7) except the 
one that has been selected. In the case of several ends the existence of such a horizon is guaranteed, 
otherwise (S,7) is assumed to have one. An area outer minimizing horizon need not be outermost, 
although the outermost one (which always exists) will have at least the same area as 5*0 (by the area 
outer minimizing property of the former). As already said, Sq need not be connected and Bray's 

theorem states Madm > \J with equality if and only if (S, 7) is isometric to the exterior region 
of the Schwarzschild manifold. 

The first step in Bray's proof is to reduce the class of metrics under consideration. To that aim. 
Bray uses an interesting result due to Schoen and Yau [142] (c.f. Theorem 7 in [30]) which states that 
given an asymptotically euclidean, three-dimensional Riemannian manifold (E, 7) (see Proposition 4.1 
in (144| for a similar statement in any dimension) with -R(7) > 0, and any positive number e, there 
always exists an asymptotically euclidean metric 7e on E, with i?(7e) > everywhere and identically 
zero outside a compact set K , which is conformally flat outside i.e. 7^ = on = Yi \ K , 
where : S°° — > R is a positive function approaching a positive constant at infinity (on each one of 
the asymptotic ends). Moreover, the following conditions are fulfilled. 



\Madm{i) - MADM{lt)\ < 



jiX,X) 



- 1 



< e, VA: e r^E and Vx G E, (38) 



i.e. the total mass changes within e and the metric itself changes pointwise at most by e (in the sense 
that unit vectors change their norm at most by e). The condition that 7^ has vanishing curvature 
scalar outside K implies that Se is harmonic (with respect to the flat metric). Therefore, the metric 7^ 
is called harmonically flat. This perturbation result has been strengthened in Bray's Ph.D. thesis [30] 
to show that can in fact be taken to be spherically symmetric, in which case 7^ is the Schwarzschild 
metric outside K. Such a metric is called Schwarzschild at infinity. This stronger version was needed 
to apply the isoperimetric methods described in Sect. 15.21 

Schoen and Yau's perturbation result implies that, in order to prove the Riemannian Penrose in- 
equality, it is sufficient to consider harmonically flat metrics. Indeed, if there existed a counterexample 
to the Penrose inequality, i.e. a metric with total ADM mass strictly smaller than the Geroch mass 
of its outermost minimal surface, then we would be able to find a harmonically flat metric 7^ with 
ADM mass within e of the original one. Since the area of the outermost minimal surface also depends 
continuously on e, it is clear that e can be arranged so that the Penrose inequality would also be 
violated for the harmonically flat metric 7c. 

For this class of manifolds. Bray defines a flow of metrics by a conformal rescaling 74 = uf-f which 
depends on a parameter t>0. The function ut is defined via an elliptic equation for vt = ^ together 
with the initial value uq — 1 (so that 70 = 7). Vt is defined by solving the Dirichlet problem 



AjVt = on Ef, 

Vf = on E \ Ef , 

vt\s, = 0, ^ 

Vt ^ — e * at 00 



(39) 



where St is the outermost minimal area surface enclosing Sq with respect to the metric 74 and E^ is 
the exterior of St. Since St is defined using the metric one is trying to construct, it is not obvious a 
priori that such problem admits a solution. Bray uses a discretization argument whereby ut{x) as a 
function of t (i.e. for fixed value of x) is discretized as a continuous piecewise linear function ul{x) 
(in t), where the jumps in the derivatives occur at fixed intervals of length e. The slope of = 
on the fc-th interval, i.e. on i S [{k — l)e, fee), fc G N is defined inductively in k. At fc = 1, is 
zero inside Sq and solves the Laplace equation A^^w^ — with boundary data = on Sq and 
—^ —1 at infinity. Sf is then defined as the outermost minimal area enclosure of Sq with respect 



to the metric 70 (the function ul^^{x) is known by continuity). Notice that since Sq was not 

assumed to be the outermost minimal surface, Si may be far away from Sq even for arbitrarily small 
e. Assume that the construction has been carried out for k G N, then the slope vl on the interval 
t € [ke, (k + l)e) is defined as being zero inside Sk and solving the Laplace equation A-ygW^ = which 
approaches the constant — (1 — e)*"' at infinity and vanishes on Sk- The surface Sk+i is defined as the 
outermost minimal area enclosure of Sk with respect to the metric (uj^^j.^j^-)^)*7o (again, the function 

'"t=(fe+i)e(2^) is known by continuity). 

Thus, while the metric is changed continuously in i, the boundaries S^, change only at discrete 
values. Existence of u'l follows from the inductive construction. The conformal flow of metrics depends 
on e. It is an important ingredient of Bray's proof to show that a suitable limit e — s- exists, thus 
defining the flow of metrics in ((39|) . The limit is performed in two ways, first the conformal factors 
are seen to converge when e — > to a locally Lipschitz function ut- This already defines the metric 
7t and the surfaces St as before. Extending the discrete collection of S^ to a one-parameter family 
via S'^{t) = SI for t E [ke, {k + l)e) the convergence of these surfaces when e ^ (for fixed t) is then 
considered. It turns out that the limit may depend on the subsequence 0. Bray denotes the 

collection of all limit surfaces as {Sa{t)}, where a identifies the element in the collection. Convergence 
is seen to be in the Hausdorff distance sense and each element in {Sa} is proven to be a smooth 
surface. The surfaces St and {Sa{t)} need not always coincide. However, they are closely related: 
for any t2 > ti > 0, St2 encloses Sa-i{ti) (for any ai) and £'02(^2) encloses St-^ (for any a2)- Thus, 
they must coincide for any value where St depends continuously on t. Bray then proves that St is 
continuous except for t in a countable set J (which may be empty). Furthermore, St is everywhere 
upper-semicontinuous and the right limit always encloses the left limit. Moreover, for t2 > ti, St2 
encloses St-^ ■ Thus, the flow of surfaces St is everywhere outward and may jump, also outwards, at a 
countable number of places. 

In addition, and this is crucial for the proof of the Penrose inequality, the area of St is constant for 
all t, even at the jumps. Outside the jumps, this constancy can be understood heuristically because 
St is a closed minimal surface and hence its area does not change to first order with respect to any 
variation, in particular the variation which transforms St into St+s {S small) in the metric 74. The 
area also changes because the metric depends on t. However, since ?j( = on St by construction, 
the metric "ft+s coincides with 74 on St to first order and therefore the area of St with the metric 74 
coincides to first order with its area with respect to "ft+s- This implies that \St\ is constant. At the 
jumps, the statement requires a careful analysis of the various limiting procedures involved. 

The second crucial ingredient in Bray's proof is the fact that the total mass m{t) of the conformally 
rescaled metric 74 is a non-increasing function of t. The proof relies on two facts. First of all, the 
definition of Vt in p9|) seems to depend on 7 and on t. However, as a simple consequence of how the 
Laplacian changes under conformal rescalings, it follows that Vt depends only on 74 (i.e. Vt solves a 
suitable Dirichlet problem with respect to the metric 74, with no reference to the original metric 70). 
Therefore, proving that < is equivalent to showing ^ | < because there is nothing that 

distinguishes 70 from any other metric 74 in the flow. For t = 0, the function vq restricted to Eq (the 
exterior of Sq in S) is just minus the Green function of 70, defined as 



The total mass for a harmonically flat metric u*S, can be computed directly from the behaviour of 
u at infinity. Namely, ii u = a + b/{2r) + 0{l/r^) then the ADM mass is m = ab. The metric 74 



is of the form 74 = (1 - tG + 0{t^)y-fo = {1 - tG + 0{t^))^U'^5, with U = 1- ^^^^ + 0(r-2), 



where Madm is the ADM mass of 79. Expanding the Green function G in the asymptotic region as 
G = 1 — c/(2r) + 0(l/r^), it follows from a direct calculation that m(i) ~ MAi:,M+t{c—2MA£iM)-\-o{t). 
Thus, the mass will not increase provided c < 2Madm- This turns out to be a general property of the 
Green function G of any asymptotically fiat metric on a manifold with compact minimal boundary 
Sq. This result, also due to Bray [22], has independent interest and has been already used in another 
context by P. Miao [117j to characterize the Schwarzschild initial data as the only static initial data 
with a non-empty minimal boundary. 

The proof of c < 2Madm relies on the positive mass theorem and uses a technique first introduced 




by Bunting and Masood-ul-Alam [53] to prove uniqueness of the Schwarzschild black hole. The idea is 
to double (So,7o) across its boundary Sq (i.e. take two copies and identify the boundaries) to define a 
new manifold E. Define also a function <i>(a;) = on one of the copies and ^{x) = on the 

other copy. This defines a function which is harmonic away from Sq (which remains a minimal surface) 
and which approaches 1 on one asymptotic end and zero on the other asymptotic end. This function 
is also everywhere. Consider the conformally rescaled metric 7 = <I>'*7o on E. The asymptotic 
behaviour of G near the infinity where it vanishes shows that (E, 7) admits a one-point compactification 
there, thus defining a complete asymptotically euclidean Riemannian manifold. Furthermore, the 
behaviour of the curvature scalar under conformal rescalings and the fact that $ is harmonic implies 
that 70 has non-negative curvature scalar away from 5*0. If (E,7) were smooth, the positive mass 
theorem could be invoked to conclude that the mass of 7 is non-negative, and zero if and only if 
(E,7) is Euclidean space. However, the metric 7 is not smooth across Sq and Bray needs to use an 
approximation argument to reach the same conclusion. Consequently, the ADM mass of (E,7) is 
non-negative. The mass of this manifold is just Madm — c/2, from which c < 2Madm follows. 

Using related techniques, Miao |116j has extended these results and has proven that the positive 
mass theorem holds for complete, asymptotically euclidean Riemannian manifolds {Q, 7) of dimension 
n < 7 and non-negative curvature scalar, with a metric which is allowed to be non-smooth across a 
compact codimension one hypersurface 5* provided (i) this hypersurface separates fl into two pieces 
ri^, with f2+ containing the chosen asymptotically euclidean end, (ii) the induced metric on S from 
both sides coincides and (iii) the mean curvature with respect to the outer direction satisfies the 
inequality 

P+{S)>p-{S), (40) 

where ^^(5*) is computed with the geometry of {n^,j\Q±). In the case considered by Bray, the 
metrics are 7± = [(1 ± G)/2]^7o. Under a conformal rescaling 7 = u'^j, the mean curvature p of a 
hypersurface changes according to p = (p + 4n(M))/u^, where n is the unit normal along which p is 
calculated. Applying this to the metrics 7± and using that Sq is minimal, it follows p^{Sq) ~ p'^{So) 
and the positive mass theorem proven by Bray is recovered. It is worth mentioning that a related 
positive mass theorem has been proven using spinor techniques by Shi and Tam 149^ on spin manifolds 
in any dimension (this entails no restriction in three dimensions, provided the manifold is orientable) 
for metrics which are smooth on ri+ and on f2^, only Lipschitz across S and such that equality in 
(jiOll holds. 

Returning to Bray's proof of the Penrose inequality, once the area of the horizons St is known to 
be constant and that the mass does not increase, the final step is to show that the metric outside 
St approaches the Schwarzschild metric in a suitable sense. The surfaces St flow outwards in E. In 
fact. Bray proves that, for sufficiently large t, St encloses any bounded set of Eq. Thus, the shrinking 
manifold Et behaves in such a way that any given point p G Eq is left out for a sufhciently large t. It 
may seem, therefore, that this manifold disappears in the limit. Here is where the asymptotic condition 
ut at infinity comes into play. For very large t, the asymptotic value of the conformal factor 

ut also tends to zero. Thus, vectors that were unit at large distances in the original metric 70 have 
increasingly smaller lengths in the metric 74. However, this metric is also asymptotically euclidean, so 
for a suitable choice of coordinates, the metric must approach the flat metric in Cartesian coordinates. 
In these coordinates, the region Et shrinks more slowly (or it may even expand). In order to see why 
is this so, consider the metric defined by 

on the domain r > i?o(i) with Ro{t) 00 as t 00. The domain is shrinking in the coordinates 
{r, 9, (/)}. However, (|4T|) is the Euclidean metric for any t. The transformation to the standard spherical 
coordinates {r',9,(j)} is given by r' — e~^*r. In the new coordinates, the domain is r' > e^^*i?o(^) 
which is not shrinking if Ro{t) grows at most as e^*. This is the type of behaviour that occurs to the 
conformal flow of metrics introduced by Bray. This can be seen explicitly by considering the flow in 
the particular case when the starting metric is the Schwarzschild manifold (|29p . Rotational symmetry 
allows one to integrate the equations easily. The flow of metrics jt in the exterior domain r > to/2 



are, for t > 0, [12]. 



_ ] ^ {dr^ + r^d9^ + sin^ Odcj)^) f <r< f e^* interior domain, 

^* ^ I (^e* + e-*)'' {dr^ + r^d9^ + sin^ Odcj)'^) r > ^e^* exterior domain. ^ ' 

The exterior domain indeed shrinks in these coordinates. However, performing the coordinate trans- 
formation (or the difFeomorphism if an active point of view is preferred) r' = re~^*, with 9 and (j> 
unchanged, the exterior region is transformed back into the original exterior part of the Schwarzschild 
manifold. The metric in the region already swapped by the surfaces {St} is cylinder-like in the 
sense that each surface r — const, has the same area. This is a general property of Bray's flow, since 
each of the surfaces St has the same area as the starting horizon. The behaviour in the exterior part 
is also general. The idea is that since the surface St grows unboundedly, eventually it swallows all the 
interior geometric features of (So,7o) and leaves only the asymptotic geometry. Since the mass does 
not increase along the flow of metrics, its limit exists and is non-negative due to the positive mass 
theorem. Bray proves that the limit mass is, in fact, positive and hence the geometry in the exterior 
approximates Schwarzschild with increasing accuracy. The asymptotic behaviour of the conformal fac- 
tor ut —^ e~* has the effect that this exterior remains unbounded, in the sense that after performing 
the appropriate diffeomorphism, all the surfaces St stay within a bounded set. More precisely. Bray 
shows that for each e > there exists a T so that for t > T there exists a diffeomorphism $t between 
(St,7t) and the exterior region r > mi?/2 of a fixed Schwarzschild manifold (Ssch)7sch) of mass mp, 
such that both metrics are e-close to each other (in the sense that the length of any unit vector of the 
metric 7* has length in the interval (1 — e, 1 -I- e) for the metric $j (7sch)). Moreover, ■m{t) and itlf 
also differ at most by e. This concludes the proof of the inequality 



Madm>\I^-^ (43) 
V Idtt 

for the original metric. As Bray points out, no property from (I],7o) inside Sq is used in the proof, 
therefore establishing (I43p also for manifolds with boundary Sq, provided this is an area outer mini- 
mizing horizon. 

While Huisken and Ilmanen's proof is strongly dependent on the dimension of the manifold. Bray's 
approach can be generalized to any dimension n < 7, as proven by Bray and Lee in |27J. The limitation 
n < 7 comes basically from two facts. First, regularity of minimal surfaces holds only for dimensions 
n < 7, and the outermost minimal surface enclosure of Sq, which is a crucial ingredient in the method, 
needs to be smooth for the argument to go through. Secondly, the positive mass theorem for general 
asymptotically euclidean manifolds (not spin) is only known to hold so far for dimensions n < 7 (the 
underlying reason is again regularity of minimal surfaces, as Schoen and Yau's proof uses this type 
of surfaces in a fundamental way). The statement of the Riemannian Penrose inequality in higher 
dimensions is, of course, different to (|43p . since the physical dimension of mass is not length in higher 
dimensions. The proper statement in (space) dimension n is 

Madai > i (^) (44) 

where is the area of an (n — l)-dimensional sphere of unit radius. The rigidity statement says 

that equality in (j44p occurs only for the higher dimensional Schwarzschild metric (first discussed by 
Tangherlini [m]) 

7^- = (l + ^^)""'5^.^-''^^^- (45) 

As discussed in [27] several of the steps in Bray's proof extend easily to higher dimensions: the 
definition and existence of the conformal flow, the proof that the area of St remains constant under 
the flow, and the fact that the ADM mass m(t) does not increase. However, the arguments required 
to show that St eventually encloses all bounded sets, and that the flow of metrics converges to a 
Schwarzschild metric after a suitable i-dependent diffeomorphism, do in fact depend strongly on the 
dimension. This is because the original argument is based on the Gauss-Bonnet theorem and a Harnack 



type inequality which is vahd only in three-dimensions. Moreover, the convergence to Schwarzschild 
requires a refinement of the rigidity part of the positive mass theorem, which roughly speaking gives 
a quantitative description of how the metric approaches the flat metric when the mass tends to zero. 
In Bray [22], this is proven using spinor techniques. This result is extended by Lee [llOj to arbitrary 
higher dimensional manifolds for which the positive mass theorem holds (in particular for dimensions 
n < 7). 

The main technical work to establish ([44]) in [27] is therefore devoted to prove that the conformal 
flow converges in a suitable sense to the metric (|45p. Besides being more general, the argument 
presented in |27| in fact simplifies and streamlines some of the original steps in the three-dimensional 
case. 



6 Penrose inequality for asymptotically hyperbolic Rieman- 
nian manifolds 

Despite its difficulty, the Riemannian Penrose inequality is considerably simpler than the general 
Penrose inequality because all the complications arising from the second fundamental form disappear. 
The time-symmetric case, however, is not the only one where this kind of simplification occurs. 
Another example is given by initial data sets (2,7^ , ^ij ) which are umbilical, i.e. when the second 
fundamental form is proportional to the metric with a constant proportionality factor. Writing Aij = 
X'fij, the Hamiltonian constraint ((Sj) becomes i?(7) = — 6A^ -t- Wnp (or R{j) = —n{n — 1)A^ + XQirp if 
S is n-dimensional). If the energy condition p > holds, then -R(7) > — 6A^ 

In an asymptotically flat spacetime, slices with this property and A ^ cannot reach spacelike 
infinity because the decay p9| is obviously not satisfied. This type of initial data are usually called 
"asymptotically hyperbolic" (also "hyperboloidal" ) since the simplest example is given by the hyper- 
boloid t = vr^ + A^ in Minkowski spacetime. As a Riemannian manifold, this hypersurface is just 
the standard three-dimensional hyperbolic space of radius 1/A (i.e. R'^ endowed with the metric of 
constant negative curvature — A^). Asymptotically hyperbolic initial data sets approach null infin- 
ity provided 7^ satisfies suitable asymptotic conditions. In the context of the Penrose inequality, 
the model example consists of the spherically symmetric slices S of the Kruskal spacetime satisfying 
Aij — Xjij. It turns out that for A > 0, E is fully contained in the advanced Eddington-Finkelstein 
portion of the spacetime and satisfies the equation (in advanced coordinates (w,r, 0, 0)) 



dv _ 
dr 



These hypersurfaces approach future null infinity for r —f 00 and intersects the white hole event 
horizon on the surface Si C {r = 2m}, which has mean curvature p — 2X (because g = 2A on any 
surface of S and 0_ = —p + q vanishes on the white hole event horizon) . The induced metric on E is 



dl^ = Jll + r^dn^ (46) 



and obviously satisfies i?(7) — — 6A^. This model space arises not only as an umbilic hypersurface of 
Kruskal, but also as time-symmetric hypersurfaces in spherically symmetric solutions of the vacuum 
Einstein field equations with cosmological constant, i.e. satisfying Ein(g) = — Ag. The general space- 
time solution is given by the so-called Kottler metric [1061 (discovered independently by Weyl [165j ) 
and reads 

ds^ = -[\ ^ + 5 ITT + r^'i'^ 47 

^ ^ r 3 

This metric is often called Schwarzschild-de Sitter (A > 0) or Schwarzschild-anti de Sitter (A < 0) 
since it contains both the Schwarzschild (when A = 0) and the de Sitter metrics (when m = 0). The 
Kottler metric admits an interesting generalization (still satisfying the vacuum Einstein equations 



with cosmological constant) whereby the spheres t — const, r — const are replaced by any compact 
connected manifold without boundary and the metric is 




) 



2m Ar^ 

r 3 



+ r'^dnl 



(48) 



where k — 1,0,-1 depending, respectively, on whether the genus of Af^ is zero, one or higher. Here 
dQk stands for the two-dimensional metric of constant curvature k. The Kottler metric obviously 
corresponds to the case k — 1. For A < and to > 0, let r+ be the only positive root oi k — 2m /r — 
Ar^/3 — 0. The hypersurface {t = 0} of (|48|) restricted to the range r > r+ > is a manifold with 
boundary which admits a smooth conformal compactification (see '38J, |40] ) to a manifold with a 
boundary consisting of two copies of AI^. One of the copies corresponds to the inner boundary r = r+ 
and the other one to the "surface at infinity", usually denoted by 9°°E. 

In the case k ~ 1 the induced metric on {t = 0} is exactly (|46|) with A = — 3A^. In this context, the 
surface 5*1 = {r = 2to} plays no special role regarding the Penrose inequality (because the surface is 
outer untrapped) and instead one has to look for minimal surfaces (similarly as in any time-symmetric 
case). Assuming again m > 0, the outermost minimal surface of (|46p is given by S2 = {r = r+j. Since 
r+ < 2m we have that Si encloses 82- 

It follows from this discussion that the Penrose inequality in this context has two flavors, one for 
which the surface of interest has p — 2|A| and another one where the surface to consider is minimal. 
The Riemannian manifolds (S, 7) relevant for this setting satisfy i?(7) > — 6A^, which is a consequence 
of ([8|) both when Aij = X-fij and p > or when the spacetime satisfies Km{g) ~ — Ag 4- 8ttT provided 
A = -SA^ and the energy-momentum T satisfies the week energy condition. 

The metric 7 has to satisfy appropriate asymptotics so that a useful concept of mass can be 
defined. Both the appropriate asymptotic behaviour and the definition of mass is considerably more 
difficult than the corresponding definition in the asymptotically euclidean case. In the case when the 
boundary at infinity is a two-sphere, the definition was first given in [163j . The definition for 

other topologies at infinity (and admitting weaker asymptotic conditions) appears in [ID] (see also [35] 
for a definition of mass from a spacetime viewpoint). For the {t — 0} slice of (|48| the mass according 
to this definition is to. 

The statement of the Penrose inequality in this setting should be expected to be different when it 
involves outermost minimal surfaces or outermost surfaces with mean curvature p = 2|A|. In the case 
of boundaries with several connected components or when the topology of is not spherical, it is 
not clear how the precise statement of the inequality should read. However, when — and the 
inner boundary S is connected, a natural version of the inequality reads 



where M is the total mass, gs is the genus of S and e = 0, 1 depending, respectively, on whether S is 
outermost minimal or outermost with mean curvature p = 2|A|. In the minimal case, this inequality 
has been proposed in [S7] for the case gs — and in [3Bj for arbitrary genus (in this reference, the 
Penrose inequality in terms of the asymptotic value of the Geroch mass under smooth inverse mean 
curvature flow is also discussed for d°°T, of arbitrary genus). In [ST], an inequality for 9°°E of arbitrary 
genus also appears. However, as noted in [35], this version fails for the slice {t ~ 0} of the generalized 
Kottler metric when has at least genus three. In the case of surfaces with p = 2|A| the inequality 
(|49|) was conjectured in [163j . Support for the validity of (|49|) comes from the fact that the {t = 0} 
slice of the Kottler metric gives equality. Moreover, by choosing the value C = — 3A^, the Geroch 
mass ([TT]) evaluated on a surface S with p — 2|A|e gives precisely the right-hand side of (|49|) . This 
value of the constant C is adapted to the inequality i?(7) > — 6A^ because then the Geroch mass is 
monotonic under smooth inverse mean curvature flow |33[ I67| . This can be seen explicitly from the 
general formula p8|l . 

Thus, it is tempting to try and adapt Huisken and Ilmanen's proof to the hyperbolic case, at least 
when the boundary S is connected. However, a recent result by Neves |130j shows that this is not 
possible in general. The difficulty lies in the fact that the Geroch mass of a flow of surfaces moved by 
inverse mean curvature does not necessarily approach the total mass of the asymptotically hyperbolic 




3 



(49) 



manifold. This is proven in [130j by showing that in the manifold {r > 2m} x 5^ with metric 
there exist surfaces with Geroch mass larger than m and which can be flowed smoothly by inverse 
mean curvature all the way to infinity. Consequently, the limit of the Geroch mass under the flow 
is still larger than m due to the monotonicity of Mg- In this example the inner boundary does not 
satisfy p — 2|A|. However, the author is able to modify the construction so that the flow starts on a 
horizon {p = 2|A|) and remains smooth for all values of the parameter in such a way that the leaves 
of the flow do not become rounder (in a precise sense) at infinity. Thus, according to the author, it 
becomes impossible to compare the limit value of Geroch mass with the value of the total mass of 
the manifold, which makes the inverse mean curvature flow inconclusive for the Penrose inequality. 
Despite this failure of the inverse mean curvature method to prove the inequality (j49p , its validity is 
still open. 

Using isoperimeric surface methods [H] (see Subsect. I5.2p . the inequality (with e = 0) has 
been proven in the special case that (5],7) outside a compact set is isometric to the Kottler metric 
(|46l) outside a sphere r = const and the following two conditions are satisfied: (i) (S,7) contains a 
unique closed and connected surface Sm with p = 2|A| and (ii) the isoperimetric surfaces Sv (with 
respect to Sm) are connected and coincide with the spheres r — const, in the asymptotic region for V 
large enough. 

It is also interesting to note that the Penrose inequality in the hyperbolic setting may be a powerful 
tool to address the uniqueness problem of asymptotically hyperbolic static initial data sets [38] . More 
precisely, denoting by U the square norm of the static Killing vector and by Mg{U) the Geroch mass 
of the level sets of U , it is proven in [38] that the validity of the Penrose inequality 




implies a uniqueness theorem for the generalized Kottler metric in the case that (?°°E if of genus larger 
than one and S is connected. The reason is that in the static setting, the Geroch mass Mq{U) can be 
bounded above in terms of the mass of the generalized Kottler solution with the same surface gravity 
K (provided this satisfies the inequality < k < |A|). The area 1 5*1 is also bounded below by the area 
of the Killing horizon of the corresponding Kottler solution. Combining these inequalities with (|5D|) . 
it follows that equality in ([50)1 is in fact the only possibility. However, the only static initial data 
set which saturates ([50]) turns out to be the {t — 0} slice of the Kottler metric. Thus a proof of the 
Penrose inequality for Mg{U) would establish a uniqueness result for static initial data sets in the 
hyperbolic setting. 

7 On the general Penrose inequality 

The validity of the Penrose inequality for arbitrary initial data sets (with an arbitrary second funda- 
mental form and without the assumption of spherical symmetry) is still open. The proofs by Huisken 
and Ilmanen and Bray of the Riemannian Penrose inequality involve manifolds S with a positive 
definite metric of non-negative curvature scalar admitting an outermost minimal surface. Although 
primarily intended to cover the time-symmetric case, the proofs only require the presence of a minimal 
surface and a Riemannian metric with non-negative scalar of curvature. Thus, the method also gives 
useful results in the case of maximal hypersurfaces, iv^A — (assuming the energy density p in ^ 
to be non- negative) , provided the Riemannian manifold contains at least a bounding minimal sur- 
face. The Riemannian arguments in the previous section, however, do not settle the general Penrose 
inequality in the maximal hypersurface case for two reasons: firstly because they would give a lower 
bound for the ADM energy instead of the ADM mass and secondly because the outermost minimal 
surface does not coincide, in general, with any of the minimal area enclosures arising in any of the 
versions of the Penrose inequalities discussed in Sect. [3] (except in the time-symmetric case, of course). 

Although the general Penrose inequality remains open, several methods have been proposed to 
address it. We devote this section to discuss them. 





(50) 



7.1 Null shells of dust 



As described in the Introduction, Penrose's original setup |136j to test the validity of cosmic censorship 
consisted of a shell of matter moving inwards at the speed of light in a flat spacetime. The shell is 
assumed to have closed and connected cross sections and the matter within the shell is made of null 
dust (meaning that the particles defining the shell are massless and that all pressures vanish). In 
order to have a flat metric inside the imploding shell, all points in their interior must be causally 
disconnected (to the past) with all points on the shell. Choosing a Minkowskian time t inside the 
shell, this demands that, to the past of some t = to, the null hypersurface Af swept by the incoming 
shell of matter must be free of self-intersections. Since the matter within the shell is moving at the 
speed of light, the cross section St = jVO {t = const} can be viewed, after the natural identification of 
points in different instants of time, as the surface lying at distance to — t from St^ , where distance is 
positive to the exterior and negative to the interior (the fact that Sf^ separates Euclidean space into 
an interior and an exterior is always true as a consequence of the Jordan-Brower separation theorem 
for smooth hypersurfaces in R", see e.g. |112) ). The distance level function from a given closed surface 
in Euclidean space is smooth everywhere in its exterior if and only if the surface is convex (i.e. it has 
non-negative principal curvatures). The setup, therefore, requires that the null hypersurface has one 
cross section St^ which is convex. This property is then true for all t < to- Towards the future, Af 
will become singular at the flrst focal point of the incoming null geodesies. A spacetime singularity 
will form there. Outside the null shell, the metric is no longer flat, in particular because gravitational 
waves may be emitted by the collapsing dust. 

Penrose devised this physical process as a potential counterexample to the inequality ([!]). The 
fundamental simplification of this problem is that the inequality can be translated into an inequality 
directly in Minkowski space, as follows [1361 1155] I157( [66] : 

Let /~ be the future directed null tangent to Af normalized to satisfy l~{t) = 1, where, as above, t 
is a Minkowskian time in the interior part of the shell. Take any closed, spacelike surface embedded 
in A/" and let be its future null normal satisfying {l^ ■ l^) — —2. The energy momentum of the 
spacetime is a distribution supported on Af which reads Ta/3 — 8'rriJ.l~l'^d, where fi is the energy 
density of the shell and the Dirac S is defined with respect to the volume form da induced by the 
normal /~ to Af, i.e. Ij^da = Ty^Q/j^e^ejeg where ef = ^p- and — » x°'{y^) is a coordinate expression 
for the embedding of N (see e.g. |125) for details). The null expansion 6*+ jumps across Af . The jump 
can be determined using the Raychaudhuri equation Q. One way of doing this is by extending the 
null vector Ij^ to a geodesic null congruence and taking its divergence on each side of the shell. This 
defines the null expansion 9+ as a discontinuous function on the spacetime. A distribution can then 
be introduced as = 0^0-1-0:^(1 — 0), where is th standard Heaviside distribution (it acts on 
tests functions by integration on the domain outside the shell) and the superscript I{E) stands for 
interior (exterior) of the shell. Since 9^0 = —l^S, the derivative of along Ij^ gives 

i'ld^e+ = iid^e"^ (1 - 0) + itd^ei® + 2 [d+] s, 

where the jump [0-|_] = {9^ — 0:^)|a/'- The Raychaudhuri equation Q, which in this case is a distribu- 
tional equation, has a singular part supported on A/ only through the term — Ric(^+,Z+) = — 327r/i<5. 
The jump must therefore satisfy 9f\j:/ = 9^_^\j^ — 167r/i. 

For an arbitrary surface embedded within a null hypersurface, the expansion along the null direc- 
tion tangent to the hypersurface coincides with the null expansion of the hypersurface. This means, 
in particular, that it only depends on the point where it is calculated but not on the specific surface 
passing through that point. This has as immediate consequence that the incoming null expansion 0_ 
of the shell is continuous across the shell. On the surface Stg , the null expansion 9- coincides with 
the mean curvature of Sto as a surface of Euclidean space with respect to the inner normal and it is 
therefore non-positive (since Stg is convex) and not everywhere zero (since Stg is closed). It follows 
from the Raychaudhuri equation that 9^ < everywhere on Af. Consequently, right after the shell 
has passed, a spacelike surface S C Af is marginally outer trapped (i.e. 9f = 0) if and only if it is 
marginally future trapped. 

Assume such an S exists along the shell. We want to show that the Penrose heuristic argument 

based on cosmic censorship then implies Mb > \J where Mb is the Bondi mass (see e.g. |160) 
for its definition) of past null infinity at the cut defined by Af . Indeed, under cosmic censorship the 



singularity that necessarily forms in the future is shielded from infinity by an event horizon. Since 
S must be contained in the black hole region, the intersection TC of the event horizon with Af must 
lie completely in the causal past of S. Using the fact that the null expansion 9- is non-positive, this 

implies \Ti\ > \S\. Since the standard heuristic argument gives Mb > \J the claim above follows. 
Notice that, in this case, the inequality is expected to hold for S irrespectively of whether or not 
this surface is area outer minimizing with respect to any spacelike slice. Using now the conservation 
equation V^qT"^ = (which holds in distributional sense, see e.g. [125j )). it follows that the integral 
/ s ^^^s does not depend on the cut 5' of A/". Evaluating this integral at past null infinity gives precisely 
the Bondi mass (this can be easily shown for instance using the Hawking mass introduced in Subsect. 
2.4p . Using now that S is marginally outer trapped from the exterior, we have, after defining 0+ =0+, 



Js Js Js 



where the last step is precisely the Penrose inequality in this setting. Thus, the Penrose inequality for 
incoming shells can be rewritten as 

/ d+Vs > yi64^, (51) 
Js 

which has the remarkable property of making no reference to the exterior geometry at all. Since 
the density fi is freely specifiable, this inequality should hold for any closed spacelike surfaces S 
in Minkowski spacetime for which the null hypersurface generated by past directed and outer null 
geodesies orthogonal to S remains regular everywhere. A similar inequality can be derived in any 
spacetime dimension |66) . 

In the particular case of a surface lying on the past null cone of a point p, the surface S can 
be described by a single positive function r which measures the distance to p (after projection to a 
constant time hypersurface) . A straightforward calculation gives the outer null expansion (?_|_ and the 
inequality ([5T|) becomes 



r7s2 > J47r / r^Vs^, (52) 



where all geometric objects refer to the standard metric of unit radius on the sphere. This inequality 
already appears in [136j and a more detailed derivation can be found in [12]. Its validity (in fact 
of a stronger version) was proven by Tod |155j as a consequence of the Sobolev inequality applied 
to a suitable class of functions on R'' (see |156| for a different proof which gives an even stronger 
inequality) . 

If S lies in a spacelike hyperplane in Minkowski, then 0+ is the mean curvature p of 5 as a surface 
of Euclidean space and the inequality becomes 

/ pris > yiHSI- (53) 
Js 

As first noticed by Gibbons in his Ph.D. thesis, this inequality is exactly the Minkowski inequality 
for convex bodies in Euclidean space, see e.g. [35]. This settles the Penrose inequality when S lies 
on a constant time hyperplane The range of validity of the Minkowski inequality (|53p has been 
extended by Trudinger [159] . to cover all mean convex bodies in Euclidean space, i.e. all closed surfaces 
with p > 0. Gibbons has used this result to claim the validity of the Penrose inequality l|5ip in the 
general case. The idea of the argument was to project S orthogonally onto a constant time hyperplane. 
The projection S can be seen to have at least the same area as S, i.e. \S\ > \S\. Furthermore, by 
direct calculation, the author finds that the mean curvature p of the projected surface is non-negative 
and that JgO+r/s = JgPTls- Thus, the Penrose inequality would follow from Trudinger's version of 
the Minkowski inequality. Unfortunately, the calculation leading to p > and Jgd+ris — JgPVs 
contains an error which invalidates the argument. Instead of going into the details of the derivation, it 
is simpler to just notice that one can easily construct a surface on a null hypersurface Af in Minkowski 
spacetime which has a projection S which is not mean convex. Consider the past null cone of a point 



a,t t — 1, and consider the sphere obtained as the cross section of Af with {t = 0}. A function s on 
S'^ taking values in [0, 1) defines a surface S on the hyperplane {t — 0} simply by moving radially 
inwards each point of the sphere a distance s. It is easy to construct surfaces S which are not mean 
convex. Consider, for instance, a surface of revolution defined by 3(6*) with equatorial symmetry and 
with a neck on the equator (i.e. such that s{9) is symmetric under 9 tt — 6 and sq = s{9 = n/2) 
is a local maximum). The principal curvatures on a point on the equator are simply 1/(1 — sq) 
and 1/(1 — Sq) + Sg/(1 — so)^, where Sq is the second derivative of s{9) at the equator. Thus, if 
Sq < —2(1 — So) (i.e. the second derivative is negative and sufficiently large in absolute value on the 
equator ) then the surface is not mean convex. But S is obviously the orthogonal projection on t — 
of the surface on Af constructed by lifting each of its points a temporal amount s to the future, see 
Fig. [21 This example shows that the projection performed in [66j is not correct. The validity of the 
Penrose inequality for null shells is therefore still open (in any spacetime dimension). 



p 




Figure 3: Surface S in Euclidean space which is non-mean convex and which can be obtained by 
orthogonal projection of a surface S lying in the past null cone of a point p in Minkowski spacetime. 



Analogously as in (|52p , the inequality (jSip can be rewritten in terms of the geometry of an arbitrary 
closed and convex surface Sq in and a function s defined on Sq, as 



(PO - 2sKo) { 1 + g^BQ^^Qj_:^fE^+l{^_zlA + ^Bq^^q^ .s(2-sPo) ) ^ 

So V (1 - + s-^K^f (1 - sPo + sH^^f ' 



> y 167r (1 - sPo + s^Ko) J7So , (54) 

where and Kq^ are, respectively, the induced metric and second fundamental form of Sq respect 
to the outer normal, Pq and Kq are the mean curvature and the Gauss curvature of S'o or, in terms 
of the principal curvatures Ai, A2 of the surface, Po = Ai + A2 and Kq = A1A2. The smooth function 
s must satisfy the inequalities s < 1/Ai and s < I/A2, but is otherwise arbitrary. It defines the 
surface S (and hence S) in a similar way as before. Details of the derivation of (|54p and some of its 
consequences will be given elsewhere [115] . 



7.2 Spinor techniques on null hypersurfaces: Ludvigsen-Vickers and Berg- 
qvist approaches 

Witten's proof of the positive mass theorem is based on the properties of spinors satisfying a suitable 
elliptic equation and which approach a constant spinor at spatial infinity. The same ideas have been 
applied to prove positivity of the Bondi mass using asymptotically hyperbolic initial data sets (see 
[4T] and references therein). Another possibility to approach the positivity of the Bondi mass is to 
use null hypersurfaces. This allowed Ludvigsen and Vickers |113j to replace the elliptic equations for 
the spinor by much simpler transport equations. Similar ideas (with different transport equations) 
allowed the same authors |114j to argue that, in spacetimes satisfying the dominant energy condition, 
the inequality Mb > ^y\S\/^&^^ holds for any weakly future trapped (in particular marginally future 



trapped) surface S which has the property that one of the two nuU hypersurfaces generated by past 
directed nuU geodesies normal to S can be extended to past null infinity while remaining smooth 
everywhere. Although this is a global assumption on the spacetime, it makes no reference to the 
future evolution of the spacetime and hence it is logically independent of cosmic censorship. For 
instance, this assumption is automatically satisfied for the incoming null shell of dust discussed in the 
previous section. However, Bergqvist |21j found a gap in the proof and the range of validity of the 
argument remains, at present, open. We discuss this next. 

Bergqvist reformulated Ludvigsen and Vickers' argument so that all spinors could be completely 
dispensed of. The idea is, in some sense, analogous to the inverse mean curvature flow for the Geroch 
mass and is based on using a functional on spacelike closed surfaces 5^ constructed as follows: start 
with a closed spacelike^surface S and let I and k hej)ast directed null normals to S partially fixed by 
the normalization {I ■ k) — —2. Assume also that I points outwards of S in the sense that the null 
geodesies starting on S with tangent vector I extend to infinite values of the affine parameter and 
intersect J^^. Let ^ be the affine parameter of this geodesic with /i = /ip on S" where /iq is a constant 
to be chosen later and assume that the surfaces = const} are smooth for all /i > /ip. The collection 
of all these surfaces defines a null hypersurface M (the "outer past null cone" of S) which is assumed 
to intersect on a smooth cut. The Bergqvist mass is defined on each leave {5^} of this foliation 
as 



provided the dominant energy condition is satisfied. Assuming that can be chosen so that asymp- 
totically 9k = —2/^ + 0{l/fj.^), i.e. without a term in then the Bergqvist mass can be seen to ap- 
proach the Bondi mass when /i — > oo. Assuming now that the initial surface is a marginally trapped 
surface, if the value of Mb could be somehow related to the area of S, a Penrose inequality would fol- 
low. Bergqvist 21] shows that this can indeed be done but only under the assumption that the induced 

S ^ ... . _9 S S 

metric y^*^ of the surface becomes round at infinity in the sense that lini^^oo M — > 9ab' 
where is the standard unit metric on the sphere. At present, it is not clear how restrictive is this 
requirement. The idea however remains interesting and deserves further investigation. 

7.3 Uniformly expanding flows 

The success of proving the Riemannian Penrose inequality in the connected horizon case using the 
monotonicity of the Geroch mass under the (weak) inverse mean curvature flow suggests a possible 
scenario for approaching the Penrose inequality in the general case. As discussed in Subsect. 12. 4i 
the Hawking mass ([T3|) is a functional on surfaces which coincides with the Geroch mass in the 
time symmetric context. Moreover, the Hawking mass (with C = 0) takes the value \/\S\/ (IGtt) on 
any topological 2-sphere which is either a marginally outer trapped surface, a past marginally outer 
trapped surface, or a generalized apparent horizon (all of which have null mean curvature vector). 
Moreover, as discussed in Subsect. 12. 4[ under suitable spacetime variations of a given surface, the 
Hawking mass is monotonically increasing [28j . Since its numerical value for large coordinate spheres 
in the asymptotically euclidean region approaches the ADM energy, a flow which interpolates between 
the horizon and infinity and falls into any of the four monotonicity cases discussed in Sect ()2.4|) would 
imply the general Penrose inequality between the ADM energy and the area of the horizon. Among the 
four cases, the closest one to the Riemannian inverse mean curvature flow is the so-called "uniformly 
expanding flow" , defined by (fTS)) . 




where dk is the null expansion along k. Using the expressions ^ and ([7]), it is easy to obtain 




where the Gauss-Bonnet theorem R{h)ris = 47rx(5') has been used. We therefore get 




With these monotonicity properties, the situation regarding the general Penrose inequahty can be 
compared to the status of the Riemannian Penrose inequahty after Geroch's heuristic argument. It is 
conceivable that the uniformly expanding flows might be useful for the proof of the general Penrose 
inequality. However, many issues remain open. For instance, in a spacetime formulation, the jumps 
that occurred in the Riemannian setting will remain, but it is unclear how and when the jumps should 
take place, even from a purely heuristic point of view. Moreover, the inverse mean curvature flow is a 
parabolic equation in the Riemannian setting, so that local existence is guaranteed, but the situation 
is quite different for the uniformly expanding flows. 

Local existence in this case has been proven only for null flows, i.e. |c| = a. A null flow will 
obviously not reach spacelike infinity. Nevertheless, assuming that a sufficiently large portion of the 
spacetime is at hand, this null flow could be used to study the Penrose inequality involving the Bondi 
mass. The difficulty, however, is that the uniformly expanding null flow seems to have the tendency to 
deform the surface in such a way that the mean curvature vector becomes causal at some places, even 
if one starts with a surface with spacelike mean curvature everywhere. This behaviour is observed 
in explicit examples in Minkowski spacetime (provided the surfaces do not lie on a constant time 
hyperplane and in fact cover a sufficiently large time interval). One alternative that may still give 
interesting results is to use a flow which is null and future directed for some interval and then continue 
with a null past directed uniformly expanding flow, then with a future directed null flow, and so on. 
This flow can be constructed in a tubular neighbourhood of the given initial data set, and therefore 
does not need strong global assumptions on the spacetime. Moreover, it is conceivable that this broken 
flow can approach spacelike inflnity. However, it is not clear which criterion should be used to stop 
the future flow and continue with the past one (and viceversa). Moreover, the construction must be 
such that the mean curvature vector remains spacelike everywhere in order to ensure monotonicity of 
the Hawking mass. 

Regarding the non-null case, |c| < a, the flow equations form a so-called forward-backward 
parabolic system, for which no local existence theory is known (this was noted by Huisken and Ilmanen 
[55] in the case c = and was extended to arbitrary c in [IH])- This is, of course, a major difficulty and 
addressing it would require a much better understanding of this type of partial differential equations. 
In Huisken and Ilmanen's work, a fundamental part of the analysis dealt with the level set formu- 
lation, which gives a degenerate elliptic equation. More specifically, the existence of a variational 
formulation turned out to be a fundamental ingredient to solve the problem of existence and to study 
the jumps. Remarkably, the uniformly expanding flows also admit a variational formulation [28]. The 
new basic ingredient is that the field to be varied is not just the level set function v (see (I36p ) as in 
the Riemannian case, but also the spacelike hypersurface E where this function is defined. Whether 
this variational formulation can give useful hints on how to define the jumps remains an open and 
difficult problem. 

It should be remarked in this context that the simplest spacelike uniformly expanding flow corre- 
sponds to c = 0, i.e. ^ — jJ^-^H [a — 1 can be chosen without loss of generality). The flow vector 

is therefore the inverse mean curvature vector. In terms of initial data information (S,7y, Ay) this 
corresponds to the case q = Q (this condition has been often termed "polar gauge" in the literature). 
Assuming this gauge condition and a second restriction which turns out to coincide with the inverse 
mean curvature flow condition pe^ = const. Jezierski was able to prove [99j llOOj the Penrose inequal- 
ity for small (but non-linear) electrovacuum perturbations of the Reisner- Nordstrom time-symmetric 
initial data outside the black hole event horizon. His method was based on writing the Hamiltonian 
constraint as a divergence term plus a non-positive reminder. The Penrose inequality is established 
by integrating this equation between the event horizon and infinity after using the Gauss theorem 
to transform the divergence into a surface integral at infinity and on the horizon (the former gives 
the ADM energy and the latter the area term in the inequality). Although Jezierski's argument does 
not use the monotonicity of the Hawking mass, the calculation does in fact correspond to the gen- 
eral identity ([TB)) specialized to the case at hand. With hindsight, one can therefore conclude that 
this monotonicity property of the Hawking mass is the underlying reason why the argument works. 
Jezierski supports the plausibility of the gauge conditions q = and pe^ = 1 by studying linear 
perturbations of Reissner-Nordstrom, where he finds that the two equations decouple, one giving a 
parabolic equation that needs to be integrated radially outwards and another one also parabolic but 
which needs to be integrated radially inwards. This is, of course, a manifestation of the forward- 



backward parabolic nature of the full system. A similar existence result for linear, axial (i.e. odd) 
perturbations of maximal slices of the Schwarzschild spacetime has been obtained in [138j . 

Another observation worth mentioning regarding the inverse mean curvature vector flow (i.e. c = 0) 
is that the monotonicity formula turns out to be insensitive to the value of the energy flux J. If follows 
that the weak energy condition is already sufficient to ensure monotonicity of the Hawking mass in 
this context. In principle this opens up the possibility that the full Penrose inequality might be 
true for spacetimes satisfying just the weak energy condition. This is not so, however. An explicit 
counterexample for scalar field initial data has been constructed by V. Husain [91] . 



7.4 Jang equation 

Schoen and Yau's proof of the positive mass theorem proceeded in two steps. First, the purely 
Riemannian case (i.e. vanishing second fundamental form) was solved by using minimal surface 
techniques. In a second step, the general case (T,,"fij, Aij) was treated by modifying the metric 7^ 
with a transformation introduced by Jang [95], namely 

Jij = Jij + dtfdjf, (55) 

where / solves the so-called Jang equation 



This transformation has the property that the curvature scalar of 7 has a lower bound that allows 
to prove existence of a conformal factor > such that the conformally rescaled metric f^^7ij has 
vanishing scalar curvature. This metric is still asymptotically euclidean and has at most the same 
mass as the original metric. The Riemannian positive mass theorem then gives the desired result. The 
proof is however involved because the Jang equation does not admit regular solutions when the initial 
data set (T,,jij,Aij) contains marginally trapped surfaces, but the idea can nevertheless be made to 
work [143]. 

A natural question, already asked in [22J, is whether a similar idea can be applied to prove the 
general Penrose inequality. The Penrose inequality has already been proven in full generality in the 
Riemannian setting so the present status is similar to the situation of the proof of the positive mass 
theorem after its Riemannian proof. This idea has been analyzed by Malec and O Murchadha [123] . 
Their argument is based on the observation that, if the Jang equation could be used to prove the 
general Penrose inequality, it should be able to do so in the particular case of spherical symmetry. 
Restricting to spherically symmetric functions /, the Jang equation becomes a simple ODE and 
existence of regular solutions in the exterior region outside the full trapped region T^t U can 
be easily shown. The difficulty of the method is that along the process (Jang's transformation and 
subsequent conformal rescaling) the mass of the manifold should not increase and the area of the 
outermost horizon S = d{T^ U 7^) should not decrease. This is because one wants to use the 
Riemannian Penrose inequality for the final manifold and conclude that the same inequality holds 
in the original initial data. It is also clear that S must be transformed into a minimal surface after 
the conformal rescaling. The simplest situation where this can be achieved is by demanding that S 
becomes minimal already for the Jang transformed metric 7^ . This requires that the outer normal 
derivative of / diverges to either +00 or to —00 on 5*. In the first case, it follows that the metric 7 has 
a cylindrical end near S. The conformal transformation is expected to compactify this end with one 
point (this behaviour was found in [143] . and this was important in order to apply the Riemannian 
positive mass theorem). Hence, the area of S decreases in the process (it vanishes in the final manifold) 
and nothing can be concluded. In the second case (outer normal derivative of / diverging to —00) the 
conformal factor is expected to have an interior local minimum. This would imply that the outer 
normal derivative of is negative on S, and hence that this surface is not the outermost minimal 
surface in the final manifold. Thus, the Riemannian Penrose inequality applied to the final manifold 
is again inconclusive for the original one. Although these arguments are not definitive in discarding 
the Jang equation method for the Penrose inequality, they indeed show that difficulties should be 
expected for the method to work, at least when S is required to transform to a minimal surface by the 



Jang transformation. The situation where S is minimal only after the final conformal transformation 
is not considered in |123j and should be further investigated. 

Very recently, the Jang equation has been successfully applied to prove a Penrose-like inequality 
in the spirit of Herzlich's inequality discussed in Subsect. 15.11 but allowing non-vanishing second 
fundamental form. The class of initial data (S, 7^ , Aij) under consideration are such that E = K(JY,°° 
with K compact and S°° an asymptotically euclidean end. The boundary dT, is compact and consisting 
of a finite collection of future MOTS (with respect to the normal pointing towards S). Moreover, no 
weakly future or past trapped surface strictly enclosing dT, is allowed to exist in S. In other words, 

is the outermost MOTS in E and, moreover, no past weakly outer trapped boundary is allowed 
to exist in S except possibly 9S itself. Assuming also that (E,7y , A^j) satisfies the dominant energy 
condition p> \J\, Khuri has recently proven [108] that the Penrose- like inequality 



holds, where k is the number of connected components of i9I] and the scale invariant quantity a is 
defined as 



The infimum is taken with respect to functions v with approach one at infinity and zero on 91]. 
The metric 7 is constructed using the Jang transformation (|55p and / is a solution of the Jang 
equation approaching zero at infinity and blowing up to +00 on dT,. The existence of such an / 
has been established by Metzger in |128j . The asymptotic behaviour of / at infinity is such that 
(I],7) is asymptotically euclidean and EADMil) = Eadm{i)- Since / ^ 00 on dT,, the level sets 
St = {/ = T}, for T large enough, form a foliation near dT which converge to dT. The idea is 
to conformally transform 7 outside St in such a way that the conformally rescaled metric jfj — 
u^Jij is still asymptotically euclidean, with vanishing curvature scalar and such that each connected 
component ST,a of St satisfies = 4^y^^/\ST,a\'yT , where is the mean curvature of ST,a with 
respect to 7^ (and the area is also calculated with this metric). Existence of ut is established as 
a consequence of the positivity properties of -R(7) and the blowing up behavior of / at dT. The 
conformal rescaling is such that the ADM energy decreases by an amount which is at least equal to 
the right-hand side of ([57)) except for terms that vanish in the limit T ^ 00. The only remaining 
step to conclude ([57)) is that the ADM energy of 7^ is non-negative. But this is precisely the content 
of the positive mass theorem proven by Herzlich [81] for asymptotically euclidean manifolds with a 
connected boundary of spherical topology and satisfying ([^ . Khuri notices that this positive mass 
theorem still holds if the boundary has a finite number of connected components of spherical topology, 
each one satisfying the bound ([5T)) . In the case at hand, this inequality is satisfied (in fact, saturated) 
by construction and the spherical topology is guaranteed by Galloway and Schoen's results [SHI HO] on 
the topology of outermost MOTS. 

Comparing this Penrose-like inequality with the difficulties described by Malec and O Murchadha 
to use the Jang transformation to prove the general Penrose inequality in spherical symmetry, the 
main difference is that this method ultimately relies on a positive mass theorem, instead of on the 
Riemannian Penrose inequality. Thus, there is no need to obtain an outermost minimal boundary 
after the metric is modified by the Jang transformation and the subsequent conformal rescaling. 

7.5 Bray and Khuri approach 

Very recently Bray and Khuri have made an important step forward towards establishing the general 
Penrose inequality. As mentioned in Sect. [3l these authors propose to use generalized apparent hori- 
zons as the appropriate surfaces for which the Penrose inequality should hold. An important guiding 
principle that led Bray and Khuri to make this conjecture is related to the fact that, independently of 
which method for proving the inequality is used, the estimates involved must all give equality when- 
ever {T,jij,Aij) is any of the slices of the Kruskal spacetime for which equality holds. Therefore, 
a preliminary issue of importance is: for which slices of the Kruskal spacetime should equality be 




(57) 




expected? Any spacelike Cauchy slice S of the Kruskal spacetime must intersect both the black hole 
and the white hole event horizons. In Kruskal coordinates, the metric reads (c.f. |160| ) 

ds^ = e^^dudv + {dO'^ + sin^ Odcf)^) , 

where uv > —1 and r{uv) is defined by uv = — 1) (our sign convention for u is different to 

that in |160| ). The null vectors dy and — 9„ are future directed, the black hole event horizon is located 
at u = 0, the white hole event horizon at u = and the domain of outer communications hi is located 
at {w > OjU > 0}. If the boundary of = 'S OU satisfies u = everywhere, then 9S°°° is an 

area outer minimizing MOTS. Moreover, since its area is lQi:m? and there are no other weakly outer 
trapped surfaces (future or past) in E"""^, it follows that this slice satisfies any of the versions of the 
Penrose inequality (|22p . ((24|l or ([25|l . It also satisfies the inequality involving generalized apparent 
horizons provided there are no generalized trapped surfaces in 1]°°'^ (this is plausible although 
not yet proven, as far as I know). Similar things happen if w = everywhere on the boundary of 
j^Doc^ However, it neither u and v arc identically zero on 91]°°^, the situation is quite different. In 
this case the boundary 91]°°° is neither a future or past MOTS and, in most cases, it is not even 
smooth. Moreover, the intersection of S with the {u = 0} hypersurface (which is always a MOTS and 
has area 167rm^) is not area outer minimizing because its mean curvature p is negative whenever v < 
(i.e. on the points lying in the white hole event horizon with respect to the second asymptotically flat 
spacetime region). Thus, its minimal area enclosure has strictly less area. Consequently, the version 
([22| of the Penrose inequality holds but not with equality (this statement assumes that {u = 0} ft S 
is the outermost MOTS of the slice, which is again plausible but not yet proven, as far as I know). 
It follows that not all slices of the Kruskal spacetime are automatically equality cases, at least for 
the version . Obviously, the more slices of Kruskal satisfy equality, the sharper is the version of 
the Penrose inequality, in the sense of being capable of identifying the Kruskal spacetime in a larger 
number of cases. A version that gives equality for any slice of the Kruskal spacetime is (pS)) . even 
when the boundary of I]°°° is non-smooth. Although, as already mentioned, no counterexample of 
this version has been found, its validity would however come as a surprise because the minimal area 
enclosure of i9(7jt U7^) may a priori have much smaller area than itself. The alternative put forward 
by Bray and Khuri involves generalized trapped surfaces. This has the advantage that, as soon as 
E°°° has a smooth boundary, this is a generalized trapped surface (in fact, a generalized apparent 
horizon). Eichmair's result, see Subsect. 12. 3( implies that an outermost generalized apparent horizon 
must exist on S. It is highly plausible that d'S°'~"^ is in fact the outermost apparent horizon in this 
case. Since its area is 167rm^, any slice with smooth 9S°°° would belong to the equality case of the 
Penrose inequality (|26p . If the boundary 5E°°° is not smooth, then it cannot be its own minimal 
area enclosure, and hence it cannot give equality neither in ((22)) nor in ((24)) . 

These considerations led Bray and Khuri to conjecture the following version of the Penrose in- 
equality (recall that our definition of asymptotically euclidean requires in particular that is 
complete and that an initial data set is called Schwarzschild at infinity if outside a compact set, 
the metric 7 is exactly Schwarzchild, c.f. the discussion after ((55)) '). 

Conjecture 1 (Bray and Khuri [29] ) Suppose that the initial data set (E,7ij, Ay ) is asymptoti- 
cally euclidean and Schwarzschild at infinity, with total mass Madm (in a chosen end) and satisfying 
the dominant energy condition p > \J\. Let S be a closed surface which bounds an open set Q con- 
taining all the asymptotically euclidean ends except the chosen one. Assume that S is a generalized 
trapped surface (with respect to the normal pointing towards the chosen end). Then 

where S'min = c^^min *s the minimal area enclosure of S (i.e. fl C fimin '"^'^ •S'min has least area among 
all surfaces with this property). Furthermore, equality occurs if and only if (S \ Q^i^jij, Aij) is the 
induced data of an embedding 0/ S \ flmin into the Kruskal spacetime such that S^i^ is mapped to a 
generalized apparent horizon. 

The use of generalized apparent horizons is indeed a completely new idea for the Penrose inequality. 
This version is not supported by Penrose's heuristic argument of gravitational collapse because it is 



not true that all generalized apparent horizons lie inside the event horizon in a black hole spacetime. 
In fact, little is known in general about this type of surfaces in general spacetimes. An important 
question regarding the plausibility of Conjecture [T] was posed by R. Wald |162j who asked whether 
there are any generalized apparent horizons in Minkowski spacetime (the question of whether surfaces 
with causal mean curvature can exist in Minkowski spacetime was already asked in [126j in a somewhat 
different context). The existence of any generalized apparent horizon embedded in a spacelikc Cauchy 
slice of Minkowski and bounding a compact set would immediately falsify Conjecture [T] Khuri has 
been able to prove [109] that no such surfaces are present in Minkowski spacetime. In fact he proves 
much more by showing that any asymptotically euclidean initial data set satisfying the dominant 
energy condition and with compact (non-empty) boundary consisting of finitely many generalized 
trapped surfaces satisfies a strict positive mass theorem Eadm > \Padm\- The proof is based on 
Wittcn's spinorial method for the positive mass. 

The strategy that Bray and Khuri propose to address Conjecture[l]is related to the Shoen and Yau's 
reduction of the general positive mass theorem to the time-symmetric case. Recall that this was based 
on the Jang transformation ([55|) of the metric, where the function / solves the Jang equation. This 
equation is specifically tailored so that it always admits a solution for any initial data set in Minkowski 
spacetime (the solution is the height function of the slice over any constant time hyperplane). This 
is relevant because slices of Minkowski immediately give equality in the positive mass theorem. A 
fundamental observation of Bray and Khuri is that, since the equality case of the Penrose inequality 
should correspond to the Schwarzschild metric, the Jang equation should be accordingly modified so 
that it holds identically for slices in this spacetime. More generally. Bray and Khuri consider static 
spacetimes {A4,g) = (K x I],ds^ = —(fp'dt^ + 7), where </> > and 7 is a Riemannian metric. The 
idea is to consider spacelike graphs {t — f{x),x) for .t G E and derive an equation for the graph 
function which holds identically in this case and which still makes sense for an arbitrary initial data 
set {'E,jij,Aij). 

The induced metric on the graph is 

which is Riemannian provided the gradient of / is not too large with respect to 7, namely if < 1. 

However, this expression can also be used to define 7^ given an initial data set (E, 7ij, A^j) and two 
arbitrary functions </> and /. This implies that the Riemannian metric 7 of any initial data set can 
be obtained as the induced metric of a hypersurface in a suitably constructed static spacetime. With 
this point of view, the function / becomes arbitrary (with no restriction on its gradient) due to the 
identity (1 — (/)^|(i/|l)(l -|- 0^|d/|^) = 1. The inverse metric is = 7'^ — v^v^ with 

where indices are lowered and raised with the metric 7ij and its inverse. The second fundamental 
form on the graph is [29] 

Of course, this second fundamental form has nothing to do a priori with the second fundamental 
form Aij of the given initial data set. Nevertheless, (E, jijjAij) is the induced geometry of a spacelike 
hypersurface in a static spacetime. The existence of the isometry generated by dt can be used to relate 
the geometry of this slice with the geometry of the corresponding {t — 0} slice, i.e. of (^,7). After 
an involved calculation, this observation leads to the following remarkable identity for the curvature 
scalar R{j) |29j . called generalized Schoen-Yau identity^ 

i?(7) = 167r(p- J,i;') A||l-|-2|z|^- |div-(0z) 

+ iY-{A - A) [ii-{A + A) + 2Aijv'v'] -f 2v'd^ {tr-(A ~ A)) , (60) 



(58) 



where Zi = [Aij — Aij)v^ and z is obtained after raising its index with the inverse of the metric 7^^-. 
In the case (f) — 1, this identity was obtained by Schoen and Yau [143] and was used to show that 
any solution / of the Jang equation defines a metric 7^ which admits a conformal rescaling with 
non-negative curvature scalar. 

Bray and Khuri's approach aims at finding appropriate functions (/> and / so that the Penrose 
inequality for (£,7^ , ) follows as a consequence of the Riemannian Penrose inequality applied to 
the transformed data (£,7). In order to use the Riemannian Penrose inequality, it is necessary that 
-^(7) > 0. The dominant energy condition p > | J|-y together with the fact that the vector v satisfies 
\v\-y < 1 (see (|58|)) implies that the first three terms in (|60p are non-negative. The remaining terms 
have no sign in general. Motivated by the structure of (pO)) . Bray and Khuri introduce the generalized 
Jang equation tT—(A — A) = or, explicitly, 



This obviously reduces to the original Jang equation ([5^ when = 1. In the present case, however, 
this equation involves two unknowns, (j) ^-nd /. The Jang equation is known not to admit regular 
solutions when S contains a future or past MOTS. In a similar fashion. Bray and Khuri observe that 
the generalized Jang equation may blow up on surfaces satisfying \p\ = which are generalized 
trapped surfaces. This, combined with the existence of an outermost generalized apparent horizon on 
(E,7y, Ay), led to authors to study the conjecture ([T]) in the particular case that dT, is a generalized 
apparent horizon and that no further generalized trapped surfaces exist in S. The conjecture in this 
setting is 

Conjecture 2 (Bray and Khuri [29j) Suppose that the initial data set (S,7y,Ay) is asymptoti- 
cally euclidean and Schwarzschild at infinity, with total mass Majjm (in a chosen end) and satisfying 
the dominant energy condition p> \ J\. Let S be a closed surface which hounds an open set E'"' contain- 
ing all the asymptotically euclidean ends except the chosen one and that S is an outermost generalized 
trapped surface with respect to the normal pointing towards the chosen end, i.e. that E = E \ E"' 
contains no generalized trapped surfaces. Then 



Furthermore, equality occurs if and only if {T,,^ij , Aij) is the induced data of an embedding of Y, into 
the Kruskal spacetime such that 9E is mapped to a generalized apparent horizon. 

Under the conditions of this conjecture, it is plausible that ([6T|) admits regular solutions in E. The 
boundary behaviour is typically singular, as examples in the Kruskal spacetime show. The hope is 
that this singular behaviour on the boundary can be adjusted so that the surface 9E has non-positive 
mean curvature with respect to the transformed metric 7^. This is useful because the area of any 
surface never decreases under the transformation 7^ 7^^- (due to the fact that the volume form of 
7y is r7=y = (1 -I- (/)^|d/|^)rj^). Consequently, the minimal area enclosure 5min of 9E in (E,7) satisfies 



where the subindex denotes which metric is used to calculate the area and the second inequality holds 
because 9E is area outer minimizing in (E, 7). Thus, an upper bound for |S'min|7 (via the Riemannian 
Penrose inequality) implies an upper bound for |9E|-y, which is the type of information relevant for 
the full Penrose inequality. 

Showing existence of solutions of the generalized Jang equation (|6ip which satisfy this criterion 
is a fundamental open issue in this approach. Nevertheless, by construction there is an interesting 
particular case where (|6ip admits solutions, namely when the initial data set (E,7ij,74y) is in fact 
a slice of a static spacetime, i.e. when there exist functions > and / such that Aij = Aij. In 
this case not only the generalized Jang equation holds trivially, but also z — 0. Using ([50)1 . this 
implies R{j) > 0. If / decays fast enough at infinity (for instance if / is of compact support) then 




(61) 




(62) 



|>S'min|7 ^ |5'min|7 ^ 



(63) 



MADAiil) — MADAiij)- Assuming that 9S is a generalized apparent horizon, it follows that 9S is 
a minimal surface in (£,7^) provided (p and / are smooth up to the boundary and 4>\di: = (it 
should be remarked that these two conditions are quite restrictive, for instance they hold for slices of 
the Kruskal spacetime only if they intersect the black hole event horizon precisely at the bifurcation 

surface {u = v = 0}). The Ricmannian Penrose inequality then gives Madm(7) > \Smin\j/ (167r), 

and hence the Penrose inequality MADAii'y) > vI^^It/O^tt) follows from ([63]) . 

Returning to the general case, whenever the generalized Jang equation is satisfied, the curvature 
scalar of 7^ reduces to 

i?(7) = 16n{p-jy) + \\A-A\\^+2\z\?^~^dWj{cl>z). (64) 

This expression has no sign in general, so the Riemannian Penrose inequality cannot be applied 
directly. However, the generalized Jang equation involves two unknowns, so it must be supplemented 
by a second condition in order to have a determined problem. Bray and Khuri discuss two possibilities. 
The simplest one consists in putting equal to zero the last summand in (|64p. i.e. 

divy((/)z ) = 0. (65) 

The two equations (|6ip -(|65 p are called the Jang - zero divergence equations in [29j . The equation (|65p 
is third order in /. However, after substracting suitable derivatives of (|6ip it can be converted into a 
second order equation for / (with quadratic second derivatives). The resulting system is degenerate 
elliptic. Bray and Khuri conjecture that the system admits solutions with appropriate behaviour at 
infinity and such that dT, is a minimal surface with respect to 7,^ . Under this conjecture, the Penrose 
inequality as stated in ([2]) follows (modulo a technical point regarding the equality case, see Conjecture 
7 in [29]). 

The second possibility is based on the observation that, while (jFTj) does not imply R{j) > 0, the 
integrated inequality (f>R{j)rpy > follows from (IM)) under suitable decay conditions at infinity and 
restrictions on dT,. Bray and Khuri make the interesting observation that in any situation (irrespective 
of whether R{j) > or not) where the mass Madm{i) can be shown to satisfy 



MadmH) ~ y > 1^ Q{x)Ria)^^ (66) 

with some Q(x) > 0, then the prescription (f> — Q imphes Madm{j) > \/|<S'min|-/(167r) and hence 
the Penrose inequality (|62p provided MADAii"/) > Madai{j)- The question is, therefore, under which 
circumstances a general type inequality of the form (|66[) holds. As the authors point out, any such 
inequality would imply the Riemannian Penrose inequality as a particular case. It is therefore natural 
to study whether the known proofs of the Riemannian Penrose inequality are capable of establishing 
the validity of for some Q > 0. The authors show explicitly that this is indeed the case for the 
Huisken and Ilmanen method, provided the second homology class of E is trivial and 9S is connected. 
The idea is to integrate (|18p with respect to A and convert the double integral (in A and on S\) as a 
volume integral. By using the weak formulation in terms of level sets, this can be accomplished even 
along the jumps. The result is that Q = \Vu\jy/ e"|S'minl7/(167r)'^/^ where u is the weak solution of 
the inverse mean curvature flow equal to zero on dT,. Thus, existence of appropriate solutions for the 
pair of equations (pT|) and (j) = Q{x) implies the general Penrose inequality ([5^ for connected dT, 
(assuming that E has trivial second homology class). As noted by the authors, existence in this case 
looks harder than for the Jang-zero divergence system because Q vanishes identically wherever the 
inverse mean curvature flow jumps. The equation therefore implies (f> — on the jumps. However, if 
/ stays smooth, then Aij = there (see (f59|) ). But then, the generalized Jang equation (j6ip requires 
tijA = along the jumps, which is a condition on the initial data and not an equation. Thus, 
existence of classical solutions of the system (p ~ Q and ((6T|) should not be expected in general. It 
may be, however, that existence can be granted if / is allowed to be unbounded (or even undefined) 
in suitable regions. 

The other existing method to prove the Riemannian Penrose inequality is the conformal flow of 
metrics due to Bray [35]. As discussed in [23], this method is also capable of producing an inequality 



of the form (|66p . In this case Q is expected to be continuous and strictly positive outside 5111111 • On 
the other hand, the resulting equation (j) = Q \s not local, in the sense that it does not define a local 
P.D.E. at each point. Existence of the coupled system with the generalized Jang equations appears 
to be difficult in this case as well. 



8 Stronger versions of the Penrose inequality 

The Penrose inequality can be strengthened under some circumstances, in particular when suitable 
matter fields are present in the spacetime. In order to understand heuristically why this is to be 
expected, let us return to the original argument by Penrose based on cosmic censorship. Assume that 
the collapsing matter is electrically charged and that the end-state of the collapse is a black hole in 
equilibrium. In this situation, all the matter sources of the electromagnetic field are expected to lie 
within the event horizon and the black hole is therefore electrovacuum in its exterior. According to 
the black hole uniqueness result, the end-state is therefore a Kerr-Newman black hole. The area of 
any cut of the event horizon in this spacetime is given by 

1^1 = 4^ (2M2 - Q2 + 2M^AP - L^/A'P - Q^^ , (67) 



where L is the total angular momentum of the final state and Q the total electric charge. From ((67|) it 
follows immediately \S\ < ATr{M + \/ — Q'^Y which makes no reference to the angular momentum. 
Since the total electric charge of the spacetime is conserved provided no charged matter escapes 
to infinity, the Penrose heuristic argument implies that any asymptotically euclidean electrovacuum 
initial data set should satisfy the inequality 



'A„in(aT+) 1 



V 167r <^li[^MADM + ^Ml^,,-Q^y (68) 

In the time symmetric case, the electrovacuum initial data reduces to the triple (S, 7^^, where the 
electric field E satisfies div^i? = 0. The total charge is defined as 47rQ = Jg{E ■ rh)ris where S is 
homologous to any large sphere in the asymptotically euclidean end. The inequality (p5)) simplifies to 

< \ [madm + v/m1^m-Q') ■ (69) 

where Sm is the outermost minimal surface. For this inequality to make sense it is necessary that 
Madm ^\Q\ for all electrovacuum initial data. This is a strengthening of the positive mass theorem 
in the presence of electromagnetic fields and was first proven in I64| (see |63j for a generalization 
including matter and |42| for a rigorous statement). The minimum value of the right-hand side of 
((69)) is IQI/2. Thus, for horizons of small area [\Sm\ < 47rQ^) the Penrose inequality ([69)) reduces to 
the positive mass theorem Madm ^ IQI, with no reference to the area of the horizon. This implies, 
in particular, that (|69p does not admit an equality case (i.e. a rigidity statement) for horizons of 
small area. On the other hand, the equality case in the Penrose inequality \Sm\ < ISttM^^j^j has the 
interesting consequence of providing a variational characterization of the Schwarzschild metric ()29p 
as the absolute minimum of total mass among asymptotically euclidean, Riemannian manifolds of 
non- negative Ricci scalar and with an outermost minimal surface of given area 15*™!. This is similar 
to the variational characterization of Euclidean space as the absolute minimum of total mass among 
asymptotically euclidean Riemannian manifolds with i?(7) > 0. 

A natural question is whether there exists another version of the charged Riemannian Penrose 
inequality which is able to give a variational characterization (among metrics of fixed charge and 
fixed area of the outermost minimal surface) of the Reisner-Nordstrom and Papapetrou-Majumdar 
metrics, which are the only static and electrovacuum black holes (see [43] and references therein). A 
strengthening of ([69)) that has been proposed is 

I I A ^ \ 

(70) 




This inequality was first discussed and proven by Jang |96| in the case of asymptotically euclidean, 
electrovacuum initial data sets (S,7ij,i?i) with a connected outermost minimal surface Sm, provided 
the inverse mean curvature flow starting on Sm remains smooth. This last requirement is, in fact, 
unnecessary thanks to the weak formulation of the flow introduced by Huisken and Ilmanen. This 
establishes (f70| for connected Sm- 

Inequality (|70p is, however, not generally true when the outermost minimal surface is allowed to 
have several connected components. A counterexample has been found by Weinstein and Yamada 
[164j . Their basic idea was to realize that the Papapetrou-Majumdar spacetime, which represents 
a static configuration of N black holes of masses nii > and charges Qi = erui, with e = ±1, has 
the property that the total area \S\ of the event horizon violates the inequality ([70)1 . Indeed, using 
Q = Madm, it follows 



irrespectively of the value of \S\. The only way how (|70p could hold is \S\ = 4'kM^jj^,j. However, 
a simple computation gives \S\ = 47r^^(TOi)^. Since the ADM mass is Madm — X^i"^*' equality 
can only happen where there is only one black hole (i.e. when the spacetime is the extreme Reissner- 
Nordstrom black hole). For any configuration with two or more black holes, ([7(I| is violated for the 
area of the event horizon. This argument is however, not a proof that ([70]) fails to hold because the 
static initial data set (i.e. the hypersurface orthogonal to the static Killing vector) in the Papapetrou- 
Majumdar spacetime does not contain any minimal surface. Each connected component of the event 
horizon corresponds to an asymptotic cylinder. Thus, some engineering is required to construct an 
electrovacuum initial data set with a minimal surface and which violates (j70p . The method followed in 
[T64j consists in an adaptation of the gluing technique developed in [93j . More specifically, it consists 
in taking two copies of a static initial data set of the Papapetrou-Majumdar spacetime with two black 
holes of equal mass mi = m2 — m, such that one of the copies has positive charges and the other one 
negatives charges. By modifying the geometry far enough along the cylindrical ends, the two copies 
can be glued together to construct an initial data set with two asymptotically euclidean ends and a 
minimal surface with two connected components. The final step is to conformally transform the data 
so that the curvature scalar vanishes. By taking m small enough, the resulting manifold violates the 
inequality (|70p . As the authors stress, this is not a counterexample of (|69p . which is the inequality 
that follows from Penrose's heuristic argument. 

Returning to the question of whether the charged Riemannian Penrose inequality provides a vari- 
ational characterization of the electrovacuum static black holes (see also [SS] for a related discussion), 
we notice that the inequality (j69[) can be written in the following equivalent way 



The Reissner-Nordstrom black holes have event horizons (or equivalently outermost minimal surfaces 
in their static initial data) with satisfy l^ml > AttQ^ (see ([S7)) with L = 0), with equality only for 
the extreme Reisner-Nordstrom case (M — \Q\)- So, these metrics belong to case (ii) above and, in 
fact, saturate the corresponding inequality. Similarly, the event horizon of the Papapetrou-Majumdar 
spacetime has area \S\ = 47rXjmf < 47r(XjTOi)^ — 4ttQ^, with equality only if there is only one 
black hole (i.e. the metric is extreme Reissner-Nordstrom again). So, this case belongs to case (i) and 
obviously the inequality is again saturated. Moreover, the Papapetrou-Majumdar static initial data 
is the only asymptotically euclidean, electrovacuum initial data (E,7y , iJ^) satisfying M = \Q\ (see 
Theorem 1.2 in 42J as well as the related previous work [63], [154| ). Thus, the formulation (|7ip would 
in fact provide a variational characterization of all static charged black holes provided the inequality 
can be proven in case (ii) with equality only for the Reissner-Nordstrom initial data. 

In the non time-symmetric case. Gibbons conjectured |65| the inequality (|70p for connected and 
outermost future (or past) marginally outer trapped surfaces S. In the non-connected case, the 
corresponding conjecture involves the sum over each connected component of the right-hand side of 
(|70p. Although no counterexample is explicitly known, such an inequality would reduce in vacuum 





if l^ml < 47rQ2 (case (i)) 
if \Sm\> 47rQ2 (case (ii)) 



(71) 



to a stronger stateraent than the standard Penrose inequahty. As noted by Weinstein and Yamada 
[164j . it seems that an initial data representing two Schwarzschild black holes sufficiently far apart 
should violate this version of the inequality. In the particular case of spherical symmetry (where S is 
automatically connected), Gibbons conjecture has been proven by Malec and O Murchadha |120j for 
maximal initial data sets and by Hayward [73j in the general case. 

In the discussion above we have dropped the angular momentum term in (|67p and have retained 
the charge. It is natural to ask what is the situation in the reverse case, i.e. when the charge is dropped 
and the angular momentum in kept. The inequality that results is \S\ < 8nM{M + \J A'P — /M'^). 
However, contrarily to the electromagnetic case, the total angular momentum of the final end-state 
after the collapse has been completed need not coincide with the initial one since gravitational waves 
carry angular momentum and this can be radiated away. As first discussed in |58| (see also |86]). there 
is one important situation where angular momentum must be conserved along the evolution, namely 
in the axially symmetric case. Under this restriction, the Penrose heuristic argument implies [79].|52|. 

A„nn(5T+) < 8^MadM [madM + \/*^ll5M " A^lm/) • (72) 

Similarly as before, this inequality only makes sense provided one can show that any asymptotically 
euclidean and axially symmetric initial data set satisfying the dominant energy condition p > |J| 
satisfies the inequality Madm > •\/F^- Again, this is a strengthening of the positive mass theorem. 
This inequality is supported by a heuristic argument based on cosmic censorship and the conservation 
of angular momentum in the axially symmetric case |54] and its validity has been rigorously proven in 
[53j for any initial data set (S,7ij, Ay) which is vacuum, maximal {iv-^A = 0), contains one or more 
asymptotically euclidean ends as well as possibly additional asymptotically cylindrical ends (which 
correspond to degenerate horizons) and such that the outermost MOTS is connected (see also [35] for 
an extension which admits furthermore non-negative energy-density and relaxes some of the technical 
requirements in [55]). Moreover, the case of equality Madm — \/\L\ occurs if and only if the initial 
data is a slice of the extreme Kerr black hole. This provides a variational characterization of extreme 
Kerr. The inequality in the case with an outermost MOTS with several connected components remains 
still open. Numerical evidences for its validity have been recently given in [55] . 

The situation for the Penrose inequality with angular momentum is therefore similar to the charged 
case. The inequality ((72| is equivalent to 



\ Ml^^, > \L\ if |5| < Svrlil (case (i)) 

\ Ml^,, > M if |5|>87r|L| (case(ii)) ^''^ 

where \S\ = ^min ((J'Tjt ) . This version of the Penrose conjecture (for axially symmetric initial data 
sets) admits a rigidity case which states that equality in case (ii) can only occur of the initial data 
is a slice of the Kerr black hole. Again, this would provide a variational characterization of the Kerr 
metric. 



9 Some applications of the Penrose inequality 

In this section we briefly mention some recent situations where the Penrose inequality has been 
exploited to derive new results. The list is probably not exhaustive, but it gives a hint on the 
potential power of the Penrose inequality as a geometric tool for adressing, a priori, completely 
unrelated problems. 

We have already mentioned in Subsect. I5.3l that the Riemannian Penrose inequality has interesting 
consequences for the quasi-local definition of mass due to Bartnik. The Riemannian Penrose inequality 
has also allowed for a dual definition of quasi- local mass due to Bray [22] (see also [26]). Here, an 
asymptotically euclidean domain with non-negative curvature scalar and whose boundary S is area 
outer minimizing is kept fixed, and all possible "fill in"s (with non- negative curvature scalar) are 
considered. The inner mass is defined as the supremum of •\/|S'|/(167r) where S is the minimal area 
needed to enclose all the asymptotically euclidean ends, except the given one. As a consequence of 
the Riemannian Penrose inequality, the inner mass is always bounded above by the ADM mass of the 
given region. The definition can also be extended to the case of non-zero second fundamental form. 



We have also mentioned in Section [5] another apphcation of the Penrose inequahty for the unique- 
ness problem of static black holes with negative cosmological constant and topology at infinity of 
genus larger than one. 

The Riemannian Penrose inequality (in fact, its proof using inverse mean curvature flow) has been 
applied recently to obtain lower bounds of the so-called Brown- York energy for simply connected, com- 
pact, three-dimensional domains with non-negative curvature scalar and with smooth boundary 
dil of positive Gauss curvature and positive mean curvature p (with respect to the outer direction). 
The Brown- York mass is defined as 



where po is the mean curvature of d^l when this surface is embedded isometrically in M.'^ . This mass is 
proven to be non- negative in [149j . Using the inverse mean curvature flow, Shi and Tam prove j l50j 
(among other things) that the inequality Msyid^) > MG{dil) holds with equality if only if is a 
standard ball in M"^. 

Still another application of the Riemann Penrose inequality is due to J. Corvino [48] . who has 
shown that asymptotically euclidean, 3-dimensional, Riemannian manifolds with non-negative cur- 
vature scalar and small mass cannot contain any minimal surface and must be diffeomorphic to R^. 
The required "small mass" condition reads 2Madm"^/K < 1, where if is a positive upper bound 
for all the sectional curvatures of the manifold. The proof is based on the fact that any outermost 
minimal surface S must satisfy \S\ > ^ as a consequence of the Gauss-Bonnet theorem. Therefore, 
the Penrose inequality implies that no minimal surface can exist in these circumstances. 

10 Concluding remarks 

In this review I have discussed the present status of the Penrose inequality. The emphasis has been 
put on trying to describe the techniques involved in the various approaches to prove it and trying to 
place the results into the right context so that a clear picture emerges of how impressive the body of 
work in this field has already been and what are the open problems that still remain. Although I have 
tried to cover the main results in this field, some aspects have been touched upon in less detail. For 
instance, I have concentrated mostly on the four dimensional case, although several results in higher 
dimensions have been discussed at various places. Some further results in higher dimensions can be 
found in the references [92], [13], [68], [5T] . 

Numerical work has also been important for a better understanding of the Penrose inequality. 
Outermost marginally outer trapped surfaces (i.e. the boundary of the outer trapped region dT^) 
are located routinely in numerical black hole evolutions in order to track the boundaries of the black 
holes. The numerical routines to do this job are collectively termed apparent horizon finders (see |152j 
for a review) and obviously they can also be used to test the validity of the Penrose conjecture. They 
have been used to check whether the Penrose inequality is fulfilled in explicit numerical examples, as 
well as for looking for counterexamples to some of its versions. They can also serve as a test to make 
sure that the MOTS being located is, in fact, the outermost one [97]. In this review, numerical results 
have been mentioned only very tangentially. Further details can be found in [102j . |103j . [90] . [52] , 



In a 3-1-1 evolution of a spacetime the outermost MOTS generates a tube of surfaces which is 
generally smooth but may jump from time to time [3J. In the smooth part, this tube has been 
called marginally outer trapped tube '32] . If the foliation is by marginally trapped surfaces (instead 
of MOTS), the tube is usually called trapping horizon ^0] or dynamical horizon [10] (with slight 
differences in their definitions). A proper study of the evolution of these tubes, specially their late 
time behaviour, is potentially a powerful method for establishing the Penrose inequality [4j. This is 
because it is expected that, at late times, the tube approaches the event horizon of the spacetime. 
Furthermore, if the MOTS foliating the tube are in fact marginally trapped surfaces, then their area 
increases with time [10 . Since this area is believed to approach that of the event horizon of the final 
black hole that forms (and this is greater than the initial ADM mass, as usual) the Penrose inequality 
would follow from a detailed understanding of the late time evolution of the spacetime and of the 
outermost marginally outer trapped tube. This approach however, is very different in spirit to the 




[I04] . [105], [97], [p; and [153^. 



ones discussed above because understanding the late time behaviour of the tube goes a long way 
towards establishing weak cosmic censorship. Thus, in essence, the Penrose inequality would follow 
because cosmic censorship would hold. So, instead of looking at the Penrose inequality as an indirect 
test of cosmic censorship, as originally envisaged by Penrose, it would become a remarkable corollary 
of a much stronger theorem establishing weak cosmic censorship (or something very close to it). In 
any case, studying the evolution of the outermost tube is an active area of research, which combines 
physical properties, numerical simulations and rigorous geometric results. The interested reader is 
referred to [1], [7j, g], [1], [140], [8j, [98j, [69j, [166], [2], [75], [3] and references therein. 

To conclude, as I have tried to show in this review, the Penrose conjecture is a very challenging 
problem that requires techniques from geometric analysis, partial differential equations, Riemannian 
and Lorentzian geometry, as well as physical intuition. The recent advances in this field have been 
impressive and our understanding of the problem is now better than ever. Nevertheless, many open 
problems remain and their study is likely to uncover new and unexpected features in the future. 
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